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Abstract 

We develop the effective theory of large-scale structure for non-Gaussian initial conditions. The 
effective stress tensor in the dark matter equations of motion contains new operators, which 
originate from the squeezed limit of the primordial bispectrum. Parameterizing the squeezed limit 
by a scaling and an angular dependence, captures large classes of primordial non-Gaussianity. 
Within this parameterization, we classify the possible contributions to the effective theory. We 
show explicitly how all terms consistent with the symmetries arise from coarse graining the dark 
matter equations of motion and its initial conditions. We also demonstrate that the system is 
closed under renormalization and that the basis of correction terms is therefore complete. The 
relevant corrections to the matter power spectrum and bispectrum are computed numerically and 
their relative importance is discussed. 
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1 Introduction 


Large-scale structure (LSS) surveys will play an increasingly important role in probing the ini¬ 
tial conditions and subsequent evolution of our universe [1]. However, extracting primordial 
information from the observations will be challenging, and understanding the many sources of 
late-time nonlinearities will be essential for realizing the full potential of the future data. On 
sufficiently large scales, fluctuations in the dark matter density are small and therefore amenable 
to a perturbative treatment [2, 3]. In contrast, dark matter is strongly clustered on small scales 
and perturbation theory is insufficient to describe its dynamics. Moreover, gravitational non- 
linearities couple short and long modes, so even the perturbative regime isn’t immune to our 
uncertainties about the evolution of LSS on small scales. 

Failing to account for the backreaction of the short-scale nonlinearities on the long-wavelength 
universe can bias the theoretical interpretation of future observations. Recently, this problem 
has been addressed using the methods of effective field theory (EFT) [4, 5]. In this approach, 
dark matter fluctuations are separated into long and short modes. While the short modes are 
nonperturbative and can only be modelled through numerical simulations, their effects at long 
distances can be captured systematically by adding corrections to the evolution equations for 
the long-wavelength perturbations. These corrections can be organized in a double expansion 
in powers of the long-wavelength fluctuations and spatial derivatives. The allowed terms are 
constrained by symmetries, and only a finite number of terms is required to describe observations 
at a finite level of precision. 

The EFT description of the long-wavelength fluctuations is complete if (and only if) the set of 
operators correcting the fluid equations is closed under renormalization. By this we mean that, 
at a given order, loops don’t generate new operators, but only mix the existing operators. The 
renormalization of the power spectrum was treated in [5, 6] (at one-loop) and in [7] (at two-loop), 
while the one-loop renormalization of the bispectrum was presented in [8, 9] (for Gaussian initial 
conditions). Extensions to halo statistics have appeared in [10-15], redshift space distortions have 
been included in [16], and baryonic effects were discussed in [17]. Since probing primordial non- 
Gaussianity is a key motivation underlying these developments, it is essential to systematically 
develop the EET approach for these more general initial conditions. In this paper, we work out 
the extension of the effective theory of large-scale structure (EET-of-LSS) to non-Gaussian initial 
conditions. We will present explicitly the renormalization of the one-loop dark matter power 
spectrum and bispectrum for a wide class of primordial non-Gaussianities (PNG). 

The importance of an accurate treatment of nonlinear corrections to the dark matter bispec¬ 
trum is illustrated in fig. 1. Shown is the bispectrum computed in standard perturbation theory 
(SPT) for PNG of local, equilateral and quasi-single-field type (see §5.3 for more details). Dashed 
lines refer to the tree-level signal, while solid lines contain the leading-order effects of gravitational 
evolution. We see that the loop corrections become relevant on relatively large scales and tend 
to decrease the difference among the three primordial signals. Even for an idealized observable 
such as the dark matter bispectrum, linear evolution is therefore insufficient to fully exploit the 
potential of current and future data sets. In this work, we will show how nonlinear gravitational 
corrections can be computed in a systematic and self-consistent way. 
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Figure 1: Solid (dashed) lines show the full (tree-level) non-Gaussian SPT contribution at 2 : = 0 
for three representative models of primordial non-Gaussianity (see §5.3 for more details). 


The outline of the paper is as follows. In Section 2, we review the shortcomings of standard 
perturbation theory and explain how they are addressed in the EFT-of-LSS. We then define the 
types of non-Gaussian initial conditions studied in this work (see also Appendix A). In Section 3, 
we show how PNG leads to new terms in the stress tensor of the EFT-of-LSS. We derive these 
terms both from the ‘top-down’, by coarse graining the equations of motion (cf. Appendix B), 
and from the ‘bottom-up’, by constructing the most general stress tensor consistent with the 
symmetries. In Section 4, we show that our set of EFT terms is closed under renormalization. 
We explicitly derive the one-loop counterterms required by the renormalization procedure (with 
details given in Appendix G). In Section 5, we compute the renormalized one-loop dark matter 
bispectrum numerically (see also Appendix D). We provide a preliminary analysis of the shapes 
of the distinct EFT contributions. A more detailed treatment will appear in [18]. Our conclusions 
are summarized in Section 6. 


Notation and Conventions 


The most important variables used in this work are collected in Appendix E. Both conformal time 
r and the scale factor a are employed as measures of time evolution. Three-dimensional vectors 
will be denoted in boldface {x, k, etc.) or with Latin subscripts (xj, ki, etc.). The magnitude 
of vectors is defined as k = \k\ and unit vectors are written as k = k/k. We sometimes write 
the sum of n vectors as fci...n = ki + ... + k^. We will often use the following shorthand for 
three-dimensional momentum integrals 


(...) 



We denote the three-dimensional Laplacian by A = S^^didj. A prime on correlation functions, 
{■■■)', indicates that an overall momentum-conserving delta function is being dropped. Our 
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convention for the dimensionless power spectrum is 

^ F{t) ^ 

We also define a dimensionless bispectrum as 

B{kuk2,h)= B{kiMM) ■ 

We will typically use p for short-scale fluctuations and reserve k for long-wavelength modes. We 
use Xi and Xg for the long-wavelength and short-wavelength parts of a quantity X. The Gaussian 
and non-Gaussian parts of X are X^ and Finally, denotes the spin-L part of X. 

When we present numerical results, the linear power spectrum is computed with the Boltzmann 
code CAMB [19], using a flat ACDM cosmology with = 0.27, = 0.73 h = 0.70. The initial 

power spectrum of the gravitational potential is taken to be of a power law form with amplitude 
= 8.7 X 10“^° and spectral index Ug = 0.96, defined at the pivot scale /cq = 0.002 /iMpc“^. 

2 Preliminaries 

We start with a review of some basic background material. In §2.1, we summarize the key 
elements of standard perturbation theory (SPT). This mainly serves to define our conventions 
and notation. For a much more detailed exposition, we refer the reader to the classic review [2]. 

Further details may also be found in Appendix G. In §2.2, we describe the effective theory of 

large-scale structure (EFT-of-LSS), and explain how it addresses some of the shortcomings of 
SPT. Finally, in §2.3, we define our way of parameterizing non-Gaussian initial conditions. 

2.1 Standard Perturbation Theory 

On large scales, dark matter behaves as a pressureless fluid, described by its density contrast 
5 = p/p — 1 and velocity v. In the Newtonian approximation, the equations governing the 
evolution of the dark matter perturbations are the continuity equation and the Euler equation: 

+ V)(I =-(l + 5)V-u , (2.1) 

[dr -I- u • V) u = —kLv — V(/> , (2.2) 

where (f> is the gravitational potential, which satisfies the Poisson equation 

. (2.3) 

Notice that the matter density parameter is time dependent. We will sometimes find it 
convenient to work with the rescaled potential = 2<f>/{3'H‘^Qrn), so that the Poisson equation 
reduces to A<I> = 6. We will also assume that the velocity is irrotational (cf. [20]), so that it is 
fully determined by its divergence 0 = V • u. 

Por small 6 and 9, the equations of motion can be solved perturbatively as an expansion around 
the solution of the linearized equations, (5(i)(fc,r) = Ili(r)5i(fc), where Di{t) is the linear growth 


4 


function, normalized to Di{Tin) = 1 at some initial time Tin (see §2.3). The nonlinear solution 
can be written as a series in powers of 5i: 

OO 

Hk,T) = J2^in){k,T) , (2.4) 

n=l 

and similarly for 6{k,T). The n-th order solution can be written as (5(„)(fc,r) ~ D^{T)5n{k), 
where 5n is a convolution of n powers of (5i: 

5n{k) = f ■■■[ {27rf6D{k-ki,,,n)Fn{ki,---,kn)di{ki)---6i{kn), (2.5) 

where = fci + ... + fcn and Fn is a kernel function which can be computed iteratively [ 2 ]. 
Diagrammatically, this kernel function can be represented by the following vertex 

= Fn{ki,---,kn){2'Kf5D{k-ki,„n)- ( 2 . 6 ) 

Correlations of 5 can then be defined in a perturbative expansion. For example, for Gaussian 

initial conditions, the power spectrum can be written as 

Ps{k) = {5{k,T)5{-k,T)y = p^i{k) + Pim + P 22 {k\ + ••• , (2.7) 

tree one-loop 

where Pmn{k) = ((5(m)(^) t))' Tperms). Notice that in order to keep later expressions 

compact we have dropped the explicit time arguments in P 5 {k) and P^nik)- Unless stated 
otherwise, density correlators should always be understood as being evaluated at time r. We will 
also be interested in the late-time matter bispectrum which, for Gaussian initial conditions, has 
the following loop expansion 

Bs{kuk2,k3) = {6{ki,T)6{k2,T)6{ks,T)y = ^ + B 222 + gfjs + -SjS + -^114 + ••• , ( 2 - 8 ) 

one-loop 

where Bi^niki, k 2 , k^) = (((i(;)(/ci, T)( 5 (m)(fc 2 ,'r) 5 („)(fc 3 , r))' + perms). The diagrammatic repre¬ 
sentations of the one-loop contributions to the power spectrum and bispectrum are shown in 
fig. 2. (See [2] for a more complete description of the Feynman rules of SPT.) 

Each loop involves an integral over an unfixed internal momentum. Conceptually, this is 
problematic since the integrals involve high-momentum (short-wavelength) modes for which 6 > 
1. Perturbation theory isn’t applicable for these modes, so there is no guarantee that the integrals 
will converge to the right answer. Indeed, this becomes manifest for power law initial conditions 
in a pure matter universe. For certain scalings of the power law, perturbation theory would 
suggest that P 13 and P 22 are divergent. The problem is less apparent in our universe, where the 
presence of radiation suppresses the contribution from the high-momentum modes. The integrals 
are then convergent, although a finite error is still being made. If the loops are only performed 
up to a finite cutoff A, then the answer depends on the arbitrary value of this cutoff scale. The 
EFT-of-LSS hxes these shortcomings of SPT. 



5 






Figure 2: Diagrammatic representation of the one-loop power spectrum and the one-loop bis¬ 
pectrum for Gaussian initial conditions. 

2.2 EFT of Large-Scale Structure 

On small scales, dark matter fluctuations are large and cannot be computed in perturbation 
theory. The EFT-of-LSS [4, 5] provides a framework in which these nonlinear modes are removed 
consistently, while keeping track of their effects on the large-scale dynamics through additional 
parameters in the fluid equations. These additional parameters cannot be computed within the 
EFT, but should be measured in observations or simulations. We briefly sketch the main elements 
of the EFT-of-LSS, but refer the reader to the literature [4-7, 20] for a more detailed account. 

We wish to describe the evolution of fluctuations above a certain cutoff scale A“^. Fluctuations 
below are integrated out, but still affect the long-wavelength dynamics through corrections 
to the equations of motion. For example, the Euler equation for the long-wavelength velocity 
receives a contribution from an effective stress tensors 

[dr +H + v^dj) u* = —(9*(/> — -djT^^ . (2.9) 

This becomes an effective theory when we expand the stress tensor in terms of the long-wavelength 
fields {4> and Vi) and their derivatives. This expansion is constrained by symmetries: in partic¬ 
ular, the equivalence principle requires the gravitational potential to appear with at least two 
derivatives acting on it and the velocity with one derivative. At first order in fluctuations, and 
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to leading order in derivatives, we have [4, 5] 


= p 


„2 Qi^ _^QiQ 


n 


n 


( 2 . 10 ) 


where the tidal tensor d'^d^ej) is absent since it is degenerate with 5*5 when we take the divergence 
of r*-^. The parameters appearing in the stress tensor (c*, c^is and Cyis) will depend on the 
cutoff A. The part of the long-wavelength solution generated by these terms will therefore also 
be cutoff dependent. This is a good thing. The new, cutoff-dependent contributions to the 
long-wavelength solution act as counterterms, in the sense that their cutoff dependence cancels 
against the one arising from loops of the long-wavelength fields. In particular, the terms in (2.10) 
are precisely the ones required to cancel the cutoff-dependent part of the one-loop contribution 
Pi 3 (fc) in (2.7). (The cutoff dependence in P 22 is removed by a noise term; see below). In fact, at 
every order in the loop expansion, the complete solution is cutoff independent if the stress tensor 
r*-^ is dehned up to sufficiently high order in the derivative expansion. This was demonstrated 
explicitly for the one-loop power spectrum in [5, 6], for the two-loop power spectrum in [7] and 
for the one-loop bispectrum in [8, 9]. So far, the discussion of renormalization in the EFT-of-LSS 
has been restricted to Gaussian initial conditions. In this paper, we extend the renormalization 
of the EFT-of-LSS to allow for non-Gaussian initial conditions. 


2.3 Non-Gaussian Initial Conditions 

We choose to define our “initial” conditions at some time Tin after matter-radiation equality, but 
early enough that nonlinearities in the prior evolution can still be ignored. The linearly-evolved 
potential 0(i)(fc,r) and the linearly-evolved dark matter density contrast 5(i)(/c,r) can then be 
written in terms of the primordial potential <p{k): 

=T{k,T)ip{k) , (2.11) 

^ii){k,r) = <t>{i){k,T) = M{k,T)<p{k) , (2.12) 

where the transfer function T(k,T) captures the evolution of the gravitational potential in the 
radiation era. To simplify the notation, we will typically drop the time arguments, but the 
functions T{k) and M{k) should always be evaluated at time r. We dehne the power spectrum 
of the primordial potential as 

p^ik) = {mp{-k)y = ^Al , (2.13) 

where we have assumed a nearly scale-invariant power law ansatz motivated by inflation. The 
power spectrum of the linearly-evolved density contrast 5(3) (A:, r) then is 

Pii(fe) = [M{k)fP^ik) . (2.14) 

Recall that both M{k) and Pii{k) are to be evaluated at time r, while Pip{k) is time independent. 
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Primordial non-Gaussianity (PNG) leads to higher-order correlations beyond the power spec¬ 
trum. The leading diagnostic for PNG is the bispectrum 

B^{ki,k2,k3) = {(p{ki)ip{k2)ipik3)y , (2.15) 

Buiiku k2, ks) = [M{ki)M{k2)M{k3)]B^{ki, k2, ks) . (2.16) 


For perturbative PNG, the potential (p can be expanded around a Gaussian field tpg. For instance, 
at lowest order, we may write [21] 

ip{k) = ipg{k) + /nl [ k-p) [ipg{p)ipg{k - p) - Pg{p) {27rf 6 D {k)] +■■■ , (2.17) 

Jp 

where we have subtracted the power spectrum of the Gaussian field, Pg{p) = {(pg{p)(pg{—p)y, 
to ensure that {(p) = 0. Eq. (2.17) is the most general quadratic expansion satisfying the re¬ 
quirements of statistical homogeneity and isotropy. The momentum-dependent kernel function 
Pvihiki, k 2 ) parametrizes the shape of the non-Gaussianity. Substituting (2.17) into (2.15), and 
keeping only the leading term in an expansion in we get 

Bg,{ki,k 2 ,k 3 ) = 2 f^^K^^{ki,k 2 )Pg,{ki)Pg,{k 2 ) + 2 perms . (2-18) 


The permutations in (2.18) imply that the bispectrum alone does not uniquely determine the 
kernel function K^i^(ki,k 2 ) (i.e. different kernel functions can give rise to the same bispectum, 
but different trispectrum). However, the ambiguity in the choice of kernel disappears in the 
squeezed limit, q = ki/k 2 ^ 1, where we get 


.^nl(^I) ^ 2) 


g-s>0 


Bg,{ki,k2,k3) 

4/NL^V(fcl)P^(fc2) 



(2.19) 


In this limit, the kernel and its leading scaling with q are uniquely defined in terms of the 
bispectrum. As we will see, the squeezed limit will be particularly relevant for our investigation. 

Statistical homogeneity requires that A'NL(fci,fc 2 ) is only a function of ki, k 2 and fci • ^ 2 - 
It will be convenient to express the angular dependence as an expansion in terms of Legendre 
polynomials Piiki ■ ^ 2)5 and write the squeezed limit as 


KNh{ki, ^ 2 ) 


(j^-O 


L,i 



A, 

Piiki ■ ^ 2 ) 


( 2 . 20 ) 


where, by symmetry, L has to be an even integer (see Appendix A for more details). The ansatz 
(2.20) captures many physically relevant cases. (Non-Gaussianity in feature models [22] may fall 
outside of this paramererization.) We can organize the different contributions by the order L 
(the “spin”) of the Legendre polynomial and treat each scaling Aj separately: 


• Scalar contributions 

First, we consider the scalar contributions; i.e. we set aL >2 = 0. The different squeezed 
limits are then distinguished by their scaling dimensions A: 
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o For A = 0, the kernel is a momentum-independent constant, K^i^{ki,k 2 ) = 1, and 
the ansatz (2.17) corresponds to a local expression for (p in real space [23], i.e. (p{x) = 
Pg{x) + f^i^{tp^g{x) — {p1))- This is the case of local non-Gaussianity. 

o For A = 2, the squeezed limit is suppressed by two additional powers of the low- 
momentum mode ki. This is characteristic of higher-derivative interactions in single- 
field inflation, which produce equilateral non-Gaussianity [24, 25]. 

o A squeezed limit with intermediate momentum scaling, A G [0,1.5], can arise if the 
inflaton interacts with massive scalar particles during inflation. In these models of 
quasi-single-field inflation [26], nonlinear interactions of the additional scalars can 
be mediated to the inflaton sector, creating observable non-Gaussianity with a char¬ 
acteristic signature in the squeezed limit: the scaling dimensions A are functions 
of the masses of the extra particles. Coupling the inflaton to operators in a con¬ 
formal field theory [27] allows to extend the intermediate momentum scaling to the 
regime A G [0, 2]. 

• Higher-spin contributions 

Various physical mechanisms can lead to an angular dependence in the squeezed limit [28]: 

o The inflaton may couple to massive higher-spin particles. In this case, the angular 
dependence is given by the Legendre polynomial of order the spin of the particle [29]. 
At tree level, only particles with even spin contribute to the bispectrum in the squeezed 
limit. 

o At loop level, the angular dependence induced by the interaction with higher-spin par¬ 
ticles can be different from just a simple Legendre polynomial. For example, coupling 
the inflaton to a U{1) gauge field via the interaction [30] leads to a shape of 

the form (2.20) with 02 = ao/2- 

o Curvature perturbations sourced by large-scale primordial magnetic helds [31, 32] also 
can lead to non-zero oq, ai and 02 . 

o The bispectrum produced in solid inflation [33] corresponds to 02 3> oq. 

3 Coarse Graining 

In the EFT-of-LSS the short-scale fluctuations have to be integrated out both in the equations 
of motion and in the initial conditions. The presence of non-Gaussianity in the initial conditions 
yields some non-trivial features which have to be taken into account in a consistent renormaliza¬ 
tion of the loop expansion. In contrast to a usual fluid (such as air or water), for dark matter there 
isn’t a hierarchy of time scales between the evolution of short-wavelength and long-wavelength 
fluctuations. This implies that the short scales keep memory of their initial conditions. For 
Gaussian initial conditions, this isn’t very important since all scales are statistically independent. 
As a result, there are no initial correlations between the long-wavelength fluctuations of the EFT 
and the short scales which are being integrated out. However, in the presence of primordial 
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non-Gaussianity, the initial statistics of the short scales depends on the long-wavelength fluctu¬ 
ations. This initial dependence will affect the dynamics of the short modes, which in turn will 
then backreact on the evolution of the long-wavelength fluctuations. The objective of this section 
is to understand how this memory effect can be incorporated in the EFT-of-LSS. 

In §3.1, we coarse grain the dark matter equations of motion to derive the stress tensor of 
the effective theory. As we explained in §2.2, this stress tensor can be expanded in terms of 
the long-wavelength helds and their derivatives. In the presence of PNG, the coefficients of this 
expansion will be spatially modulated. In §3.2, we derive this dependence by coarse-graining the 
initial condition (2.17). This introduces new, non-dynamical fields in the EFT. We collect the 
leading-order non-Gaussian corrections in §3.3. 

3.1 Smoothing the Equations of Motion 

We will follow the notation of [4] , and write the long-wavelength fluctuations of any held X as 

Xe(x) = [X]a{x) = I d^s' Wa{\x - x'\) X{x') , (3.1) 

where W\{x) is a window function which quickly vanishes for x > A“^. In momentum space, 
this convolution becomes a simple product 

X,{k) = FA{k)X{k) , (3.2) 

where F^ik) is the Fourier transform of W\{x). Naturally, the short-scale huctuations of X are 
dehned as the complement Xg = X — Xg. 

While the continuity equation remains unchanged,^ the smoothing of the Euler equation (or, 
more generally, the Vlasov equation) generates terms which contribute to an effective stress 
tensor [4, 5, 34] 

- V^(5fc((),)2]^ + [pvlvi]^ + [pdf]^ , (3.3) 

where is the velocity dispersion of the short scales. To determine the effect of this stress 
tensor on the long-wavelength universe, we need to understand how correlates with long- 
wavelength modes. When the short modes are perturbative, this backreaction on the long modes 
can be computed explicitly in perturbation theory (see Appendix B). However, in general, the 
short scales are not perturbative and therefore cannot be computed analytically. When this is 
the case, the stress tensor should be written as the most general function of the long-wavelength 
fluctuations consistent with the symmetries of the problem. To do so, let us first notice that 
the dynamics of the short-scale fluctuations depends both on their initial conditions and on the 
long-wavelength background in which they evolve. Assuming that the short scales depend locally 
on the long-wavelength fluctuations, the most general effective stress tensor consistent with the 
equivalence principle then is 

F^{x,t) = F[(f)s{q,Tin) ■, d^d^l{Xfl{T'),T') , d^vl{Xfl{T'),T') , •••] , (3.4) 

^This is only true if we smooth the momentum density field tt = pv instead of the velocity field v [20]. 
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where is the position of the fluid element which at time r > r' is located at a;. It 

will be important that the initial short-scale fluctuations are evaluated at the Lagrangian posi¬ 
tion q = Xfi(Tin)- The ellipses refer to higher derivatives in the long-wavelength fluctuations. 
Since there is no hierarchy in the time scales of the evolution of the short and long modes, the 
stress tensor depends on the evolution (along the fluid trajectory) of the long-wavelength fluc¬ 
tuations. However, as explained in [8], when perturbation theory is valid, this expansion can 
be reorganised in an expansion which is local in time, but in general non-local in space (unless 
convective time derivatives are added to the expansion, see [35]). Working at second order in 
perturbation theory, one can show that (5(®j;(r'), r') is locally related to the late-time overden¬ 
sity 6{x,t) and the tidal tensor didj(f>{x,T).‘^ Hence, in this work, we may assume that the 
dependence of the stress tensor on the long-wavelength fluctuations is effectively local in time. 

At first order in the long-wavelength fluctuations and to leading order in derivatives, we can 
write (3.4) as 

T*-^’(£c,r) = c*-^(s,r) -h d^kiix,T)d^d’'(j)i{x,T) + ki{x,T) d^v\{x,T) + ■■■ . (3.5) 

We will sometimes suppress the indices on the coefficients in (3.5) and collectively refer to them 
as c{x,t). Notice that the coefficients c{x,t) depend on position through their dependence on 
the initial short-scale fluctuations^ 


c{x, t) = c[cl)s{q, Tin),T] . (3.6) 

Since the short scales have been removed by coarse graining, the relation in (3.6) is non-local. 
More precisely, the coefficients c{x, r) will depend on the value of cj)s over a patch of size 
centred around the Lagrangian coordinate q{x,T). We will come back to this point in the next 
section. 

When the small-scale fluctuations are replaced by their statistical ensemble averages, the stress 
tensor becomes 

+ , (3.7) 

where ■ ■)s denotes an average over many realizations of the short modes in a fixed long- 
wavelength background and Ar*-^ is a stochastic term which accounts for the statistical deviation 
from the average. We will sometimes refer to as the “viscosity” part of the stress tensor and 

to Ar*-^ as the “noise” part. We get the viscosity contribution to the stress tensor by replacing the 
coefficients c{x,t) in (3.5) by their statistical averages {c{x,t))s- For Gaussian initial conditions, 
the average over short-scale fluctuations is independent of the long-wavelength fluctuations. As 
a result, the averaged coefficients {c{x,t))s become simple cutoff-dependent parameters of the 
EFT. However, primordial non-Gaussianity gives rise to correlations between different scales and 

^ Note that this relation relies only on the equations of motion and therefore applies regardless of whether the 
initial conditions are Gaussian or not. 

^In (3.6), we have assumed that the initial short-scale fluctuations are fully determined by the short-scale 
potential (jisfq, Ti„). We could also have added a dependence on the short-scale velocity vl{q, Tin) and its dispersion. 
However, at early times, every scale is perturbative so that one can use the linear equations of motion to express 
the velocity in terms of <j)a- In particular, assuming that these scales are in their growing mode, one can show that 
the velocity is determined by the gradient of the potential, w* oc d^(j)s- 
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the averaged coefficients will depend (non-locally) on the initial long-wavelength fluctuations. 
This dependence can be determined by integrating out the short-scale fluctuations in the initial 
conditions, which we shall do next. 


3.2 Smoothing the Initial Conditions 

We apply the filtering procedure (3.1) to the primordial potential (2.17). Isolating the terms 
which contribute to the coupling between the long and short modes, we get 


^k) D /nl / K^i^ik - p,p)[<pl{k - p)ip^g{p) - Pg{p) { 27 Tf 6 D{k)] = fNi.i^j{k) , (3.8) 

(fsip) D 2/nl [KNl.ip-P,P)^gip-p)^'gip) ■ (3.9) 

Jp 

The field V’J in (3.8) does not correlate with the long-wavelength fluctuations, but it has non¬ 
trivial correlations with the noise terms of the stress tensor, Ar*-^. We will explain this in more 
detail in the next section. The right-hand side of (3.9) encodes the dependence of the short modes 
on the long-wavelength fluctuations. Writing 4>g{p,Tin) = T{p,Tin)(Pg{p) and going to position 
space, we find 


4)s{x,Tin) ^ 4)t{x,Tin) +2f^ 


k Jp 




(3.10) 


Since k p, the result in (3.10) only depends on the squeezed limit of the kernel function 
(and the corresponding bispectrum). In §2.3, we introduced the following ansatz 


K^^{k,p) 


Piik-p) l + 0{k‘^/p^) 


(3.11) 


For each scaling A in this sum, we can treat the different orders (“spins”) in the Legendre 
expansion separately: 


• Spin-0 

For L = 0, the initial short-scale fluctuations can be written as 

(t>s{x,Tin) = 4>l{x,Tin) + fNl.a^^\x)'ll;{x) , (3.12) 


where 



(3.13) 


and a^^\p) = 2ao (p/p)^ (t^g{P: Tin)- The scale p is an arbitrary momentum scale introduced 
to make ipix) dimensionless. As we shall see, the field V' will play an important role in the 
rest of this paper. It describes how the statistics of the small scales is modulated by the 
presence of the long modes. For local non-Gaussianity, we have A = 0 and therefore 
'4>{x) = (pg{x). Eq. (3.12) then reduces to a more familiar expression 



(3.14) 
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We see that local non-Gaussianity simply modulates the amplitude of the small-scale fluc¬ 
tuations. 


To compute the dependence of the coefficients (c(£c, t))s on ■0, we substitnte (3.12) into (3.6). 
At first order in an expansion in ip, we find 

Sc{x, r) 


{c{x,t))s = {c[(l)i{q,Tin),T])s + h 




(3.15) 


)q \^4’s{qjTin) 

where the derivative is evalnated for the Gaussian field configuration cpg. Furthermore, 
since the integral only has support for \q — q{x, r)| < A“^, we can pull the long-wavelength 
field ip{q) ~ P^iq) out of the integral (up to corrections that are higher order in V^/A^). 
Translation invariance then gnarantees that the remaining integral is independent of the 
position X, and eq. (3.15) can be written as 

{c{x,t))s = c(A,r) -F /NLC.i/,(A,r)T(a;,r) , (3.16) 

where we have defined 

^{x,T) = ip{q{x,T)) . (3.17) 

In Section 4, we will see why it is important that field ip is evalnated at the Lagrangian 
position q. Restoring the indices, but dropping the arguments of the coefficients and the 
fields, we get 

, (3.18) 

= [cf' + /nl + [cfl + U , (3.19) 

and similarly for the coefficient The tensor structure of (3.18) and (3.19) is fixed 

by statistical isotropy. 

Spin-2 

For L = 2, equation (3.10) becomes 

(ps{x,Tin) = (pg{x,Tin) -b /nl (®) (®) 

where 


(3.20) 




(3.21) 


and a\j {p) = 02 {pi/p)^ {PiPj)(pg{P^'''in)- In (3.21), we have defined the projection operator 
V^^{k) = k'^k^ — Note that the tensor ip'^^ is traceless and hence contains no scalar 

contribution. Substituting (3.20) into (3.6), we find 

+ , (3.22) 




(3.23) 


+ /n 


b-o 
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where we have defined r) = r)). Although we have dropped the arguments, 

the coefficients in (3.22) and (3.23) of course still depend on the cutoff A and time r. 


• Spin-4 

For A = 4, equation (3.10) becomes 

(psix,Tin) = (pl{x,Tin) + fj,i^a[fj^i{x)'lp^^’^\x) , 

where 



(3.24) 


(3.25) 


and alfj^iip) = 04 {p./p)^ {PiPjPkPi)4>g{P^ Tin)■ The projection tensor ensures that 

^ijki -g gynimetric and traceless in any two of its indices 

(^5''^k^U + 5 perms 

Let us note that, at leading order in derivatives, the coefficients with only two indices, c*-^, 
cannot depend on since we would need to contract two of the indices. As a result, the 
field only contributes to the coefficients with four indices 

W = + /nl , (3.27) 

where, as before, we have defined = 'il)^^^\q{x,T)). 


= yyy^u - ^ 


^ ^ + 2 perms^ . (3.26) 


• Higher-spin 

It should be clear that higher-spin contributions (A > 4) will be parametrized by higher- 
order symmetric traceless tensors. However, at the order we are working at, these higher- 
order tensors will need to be contracted with one or several Kronecker deltas. Hence, they 
do not contribute to the stress tensor and will therefore not be considered in this work. 


3.3 Effective Stress Tensor 

Putting the results of the previous two sections together, the stress tensor becomes a functional 
of the long-wavelength fields^ = 20/(37^^11^)}) the non-dynamical fields {'k, 

and their derivatives. In this section, we present the most general expression for the viscosity 
and noise parts of the stress tensor. 

Viscosity part 

We write the viscosity part of the stress tensor as 

(r^i)^ = d¥, T, • • • ] , (3.28) 

^To avoid clutter in the expressions, we will drop the subscripts on the long-wavelength fields from now on: 
i.e. we set {v\, 4>^} —^ {u*, I?}. 
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where the ellipses refer to terms with higher derivatives. Since the long-wavelength fluctuations 
are perturbative, we can expand the functional J- in (3.28) in powers of the fluctuations and 
their derivatives. We will focus on terms which involve single powers of the new long-wavelength 
fields since the Gaussian terms have been discussed in [5, 8 , 9, 36]. The stress 

tensor can then be written as the sum of three terms organized by spin: 





(3.29) 


We look at each of these contributions in turn. 


• Scalar contributions 


The leading spin-0 contributions to the stress tensor are 


1 

P 



+ ci^S + ^^6 

ri 


5^^ + C3 + d^v^) , 

2'rL 


(3.30) 


where the coefficients ^(r) and Ci(r) are dimensionless and time dependent. At second order 
in perturbation theory, the expression in (3.30) simplifies, since some terms are related by 
the equations of motion. In particular, using with / = dlnDi/dlna, 

we get 


= [g^ + + g2-^s^^ , 


(3.31) 


where we have defined = ci -|- |c 3 — / C 2 , 52 = C 3 — / C 4 and introduced the tidal tensor 
Sij = didj^ - A$. 


• Higher-spin contributions 

Similarly, the leading spin-2 and spin-4 contributions to the stress tensor are 


p ^ 


(3.32) 

(3.33) 


We see that a bispectrum with a squeezed limit of the form (3.11), generally yields eight additional 
parameters in the stress tensor (for each value of Aj). When the squeezed limit is isotropic, the 
number of additional parameters reduces to three. 


Noise part 


The noise term in (3.7) arises because, for any specihc realization, the short-scale modes fluctuate 
away from their ensemble averages. This stochastic contribution to the effective stress tensor can 
also be expanded in powers of the long-wavelength fields {4>, T, , 4'*-^^^}. Just like the viscosity 
term in (3.29), the non-Gaussian contributions to the noise term can be organized in terms of 
spin: 




(Ar*J)[oi + {AHip] + (Ar*^')W 


(3.34) 


We look at each of these contributions in turn. 
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• Scalar contributions 

The scalar contribution is given by 


(3.35) 


where is a random variable representing stochastic noise, i.e. contributions which are 
uncorrelated with the long-wavelength fluctuations. This noise term also correlates with the 
leading-order Gaussian noise term, D Jq . As we shall see in Section 4, the correla¬ 

tion of these two noise terms will be important in the renormalization of non-Gaussian loop 
diagrams. Furthermore, noise terms are uncorrelated on large scales, i.e. their correlation 
functions in position space are proportional to delta functions. This means that, in Fourier 
space, the correlator of with is 

{ji\k,T)4\k\T)) = J^^^\k,T){2TTf5D{k + k') , (3.36) 


p ^ 


where r) is an analytic function, which around k = 0 can be expanded as 

t) = J{A, + J{A, t) + 0{k‘^) . (3.37) 


Furthermore, as advertised earlier, the noise terms of the stress tensor, Jq^, also correlate 
with the noise term in the initial conditions, ijjj'. 


{i;j{k)J^J{k',T)y = j K^^{k-p,p){ipl{k-p)(pl{p)Jiy{k',T))' 

Jp 

= Ja{A,T)S^^ + 0{k^), (3.38) 

where a G {0,ip} and we have assumed that the integral is analytic in k. 

• Higher-spin contributions 

The leading spin-2 and spin-4 contributions to the noise part of the stress tensor are 

hAT«)Pl = jf‘<i!a , (3.39) 

tjAC’P' = , (3.40) 

where the tensors and are uncorrelated with long-wavelength fluctuations, but 

can be correlated with and Jq . 


Both the viscosity and noise parts of the stress tensor are crucial in a consistent renormalization 
of matter correlation functions. As we will see in the next section, they generate new solutions 
which have precisely the correct momentum dependence to absorb the divergences arising from 
the loop diagrams of standard perturbation theory. 
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4 Renormalization 


In this section, we describe the one-loop renormalization of the power spectrum and the bis¬ 
pectrum in the presence of primordial non-Gaussianity. We will focus on loop diagrams arising 
from non-Gaussian initial conditions, since loops coming from Gaussian initial conditions have 
already been studied in [7-9] . Moreover, we will restrict the presentation to the case of isotropic 
(spin-0) initial conditions. The generalization to PNG with non-trivial angular dependence is 
straightforward, and left as an exercise to the reader. 

In §4.1, we classify the divergences arising in the non-Gaussian loops of standard perturbation 
theory. We then show, in §4.2, that the solutions generated by the stress tensor (3.31) act as the 
appropriate counterterms. When we compute correlators of the field 'I' it will be convenient to 
expand the right-hand side of (3.17) around the Eulerian position x \ 

4'(a;,r) = 

= if^ix) -|- VV'(a;) • V4>(®, r) -!-•••, (4.1) 

where, in the last line, we have used linear perturbation theory to write the velocity n* in terms 
of the gradient of the gravitational potential 9*4. The correlation between 5 and ijj is 

PiV,(fc) = ((l(i)(fc,r)V>(-fc))' = i^^Pn(A:) , (4.2) 

where A and M{k) were defined in (3.13) and (2.12), respectively. Finally, in §4.3, we discuss 
the time dependence of the renormalized EFT parameters. 


4.1 Loops in Standard Perturbation Theory 

There is only one non-Gaussian contribution to the one-loop power spectrum 


Pi2{k) = 



/ F 2 {p,k - p)Biii{k,p,\k - p\) + perm 
J n 


Extracting the UV behavior of this integral, we find the following divergence 


(4.3) 


Pi2{k) 


p^oo 


where we have defined 


-^/nl [2a‘^{A)k^ Pi^{k) + a^{A) k‘^M{k)] + perm , 



dp gp 
27r2 (p/p)'^ 


-Pii(p) , 


dp K^^{p,-p) 
27r2 [M(p)]2 


Pfiip) • 


(4.4) 


(4.5) 

(4.6) 


Recall that M(p) and Pii(p) are evaluated at time r (or scale factor a). The coefficients <7^(A) and 
cj|,(A) are therefore also time dependent. In (4.6), we assumed that K^i^{p,—p) is independent 
of the direction of the loop momentum p. We see that the first term in (4.4) is proportional to 
Pi^{k), and will therefore be renormalized by a counterterm proportional to ip. On the other 
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hand, the second divergence in (4.4) is proportional to M{k) and is renormalized by the noise 
contribution to the initial conditions V’j- 

There are four non-Gaussian contributions to the bispectrum (see fig. 3): 


-^113 = 3-Pii(A;2) [ Fsiki + p,-p,k2)Bniiki,P,\ki + p\) + 5 


perms 


fif/] = 3Biu{ki, k2, ks) f F3{ki,p, -p)Pii{p) + 2 


perms 


-^122 = ‘^ F 2 {k 3 +p, -p)F 2 {p,k 2 -p)Bui{ki, 1^3 +p|, 1^2 -p\)Pii{p) + 2 perms , 

Jp 

-^122 = F 2 {ki,k 2 )Pii{k 2 )Pi 2 {ki) + 5 perms . 


(4.7) 

(4.8) 

(4.9) 
(4.10) 



Figure 3: Diagrammatic representation of the non-Gaussian contributions to the one-loop bis¬ 
pectrum. The diagrams of type II are renormalized by the same counterterms that renormalize 
the one-loop power spectrum. 


We note that the divergences appearing in B ^22 already arise in P 12 and are therefore renormalized 
by the same counterterms. Hence, we turn our attention to the diagrams and Bf 22 - 

• The UV limit of is 


B 


(I) 

113 


p^oo 


4 ,2 , 1 ki k 2 2 , 7 {k3-k2f 

105 ^^ ”^21 kl ^”^45 kl 


X /n 


2cT2(A)Pi^(fci) + 4(A)M(fci) 


Pii{k 2 ) + 5 perms 


(4.11) 


where (T^(A) and o'^(A) were defined in (4.5) and (4.6), respectively. Let us focus on the 
term proportional to fci • ^ 2 - It is easy to see that this contribution is renormalized by the 
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counterterm Vi/’ • V$. However, in order for this term to transform correctly under boosts, 
it has to appear in the combination 


'(/’(a:) + VV’(») • V$(a;,T) , 


(4.12) 


i.e. the field i/’ has to be evaluated at the Lagrangian position q{x, r); cf. eq. (4.1). For this 
to happen, the divergences in (4.4) and (4.11) have to have the same numerical coefficient 
(highlighted in red). It is reassuring that our computation reproduces this fact. 

• The UV limit of is 


B 


(ii) 

113 


p^oo 


61 

= -^^'(A) {kl + kl + kl) Hin(A;i, k^, k^) 


(4.13) 


where we have defined 

f^^(A) = J ■ 


(4.14) 


Since the divergence in (4.13) is the same as that appearing in P 13 and P 123 (cf. fig. 2 ), it 
will be renormalized by the same (Gaussian) counterterm. 


• The UV limit of B ^22 


B 


(I) 

122 


p^oo 


= -—f 

147 /NL 


kjkl - ^(fc2 • kaf 


2&\A)P,^{ki) + al{A) M{ki) 


+ 2 perms , 


where we have defined 

d 2 (A) 

(A) 



— p 2 („) 

27r2 p2(p/^)A^llW ’ 

^ K ^^ jp , - p ) 3 

27r2 p‘^[M{p)]‘^ 


(4.15) 


(4.16) 

(4.17) 


We notice that all terms in (4.15) are analytic in two of the external momenta. For instance, 
the first permutation is analytic in k 2 and k^. In position space, these divergences therefore 
are proportional to derivatives of delta functions. As we will explain below, these terms are 
renormalized by the noise term in the effective stress tensor (3.35). 


4.2 Renormalization in the EFT-of-LSS 

Next, we demonstrate that the stress tensor (3.31) generates new solutions which precisely cancel 
the cutoff dependence arising from the loop diagrams in SPT. Many of the computational details 
will be relegated to Appendix C. 

Including the stress tensor, the equations of motion for 6 and 9 are 


Vs5 = 


—0^9? + ( —~ 3 ) ada + 


6 = Si3 - BdaiaSa) + Te 


VgO = 


+0^9^ + ( 4 — -Hm ) ada + (2 — 30^) 


9 = da{aSp) - -VLmUSa + da[aTe) 


(4.18) 

(4.19) 
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(4.20) 


where So, and 5^ are nonlinear source terms defined in Appendix C and 


Te = -di 


P 


Substituting (3.31) and (3.35), we find that tq is the sum of a contribution from {t^^)s and one 
from the stochastic term 

re = Ty + Tn , (4.21) 


where, to second order in the fluctuations and including the Gaussian contributions (see [8, 9]), 
we have 


Ty = —<f A5 — ei A((5^) — 62 A(s^) — e^di{s'^^djS) 

- - di{5d^^)) + + g2didj{^s^^)] , (4.22) 


Tn = -di [djji^ - {5djji^)\ - didj{J\^6) - didj{f^MS^^) - . (4.23) 

The full solution {<5, 6} can then be written as a sum of three terms 

5 = jSPT Q ^ ^SPT _^qc^qJ ^ (4 24) 

where is the SPT solution (see Appendix C), while {S'^,9'^} and {6"^,9'^} are the 

solutions generated by and r^,, respectively. In (2.4), we expanded the SPT solution in powers 
of the first-order initial condition (5i. A similar expansion can be defined for the counterterms 

OO 

= , (4.25) 

n=l 

where oc (^i)”" and (x ((Ji)”"^, and equivalently for 9. In particular, the viscosity 
counterterm can be expressed as the sum of a Gaussian and a non-Gaussian contribution 

(5A(fc,o)= / •••/ (27r)^(5D(fc - fci...n) • ,fcn|a) 5(i)(fci,o) ■ • •5(i)(fcn,o) 

Jkl Jkn 

+ /nl /" ••• /" {2Trf5D{k - ki,„n) .. ,kn\a)'ijj{ki).. .6(^i){kn,a) , (4.26) 

J t/ 

where the kernel functions and have the following diagrammatical representations: 

= F^{ki,...,kn\a){2TTfSD{k-ki...n), (4.27) 


= UH^{ku...,kn\a){27rfSD{k-k,,„n) . (4.28) 
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Our goal in this section is to compute and 6'^ up to second order. Naturally, both <5'^ and 
have contributions from the Gaussian and non-Gaussian terms of the stress tensor. However, 
since the Gaussian contributions have already been computed in [ 8 , 9, 36], we will mainly focus 
our attention on the non-Gaussian contributions. Furthermore, since we are only interested in 
showing that the counterterms cancel the cutoff dependence arising from the SPT loop diagrams, 
we will, for simplicity, restrict the analysis to the Einstein-de Sitter cosmology. The extension to 
AGDM will be explained in more detail in the next section. 


Viscosity counterterms 

• The first-order counterterm is the solution to 

Vs 61,) = -d\a) A 5 ( 1 ) - g{a)U A V' • (4.29) 

The Green’s functions associated with the operator Vs is defined in Appendix C. It allows us to 
write the solution of (4.29) as 


6l,){k,a) 

= -^(a)A:^5(i)(fc,a) - 7(a)/NLfe^V’(fc) , 

(4.30) 

where we have defined^ 



C{a) 

= — f da' Gs{a,a') a'd^{a) , 

(4.31) 

7(a) 

1 = - / da' Gs{a, a') g{a') . 

J Clifi 

(4.32) 


The non-Gaussian contribution to ^l’^) cancels the loop divergence in P 12 : 


P12 +Pic ^ ^ 


-|- □—^-•-►--H 


= -/nl + 27 («)^ k^Pi^{k) = finite , 


(4.33) 


where, for the moment, we have only focused on the divergence proportional to Pi^(A:). (The 
divergence proportional to M{k) in (4.4) will be cancelled by a noise term.) Hence, the cutoff 
dependence of the one-loop power spectrum is removed if 

2 1 

7 (A,a) =- c 7 ^(A, a)-|-finite 47 g{A,a) = -C 7 ^(A, a)-|-finite , (4.34) 

21 3a 

where we have added the explicit arguments to highlight that both 7 and a are functions of the 
cutoff A and the time a. The finite piece of (4.34) will be discussed in §4.3. 


®Note that the parameter 7 of [ 8 ] corresponds to the parameter ^ in this paper. 
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( 11 ). 


The Gaussian contribution to cancels the loop divergence of 


B 


(ii) 

113 




+ C(®) (^1 + ^2 + k‘^)Bui{ki, k2, ks) 


Hence, renormalization requires 
61 

({A, a) = 0 - 2 (A, a) + finite 

DOU 


61 1 

d^(A,a) =- —-cj^(A,a)+ finite . 

ioU a 


(4.35) 


(4.36) 


Note that the parameter ^ is fixed by the renormalization of the Gaussian one-loop diagram P 13 [5] . 


The second-order counterterm 5^2) satisfies 


:( 2 ) 


^5<5(2)=4 -'^'9,(a5(2)) + r(2) , 


(4.37) 


( 2 ) 

where ^ are obtained by replacing one of the 6 (or 9) in the nonlinear source terms by their 
first-order SPT solution and the other one by the first-order counterterm '^( 1 ) ^( 1 ) = -'^( 1 ))- 

Focusing on the solution generated by the non-Gaussian terms, we hnd 

(5f’^'^(fc,a) = [ da'Gs{a,a') S^g\a') - nda/ia'S^^a')) + {a' 


42 ) 


= /nl / i7|(p,fc-p|a)V’(p)-5(i)(fc-p,a) . 

Jp 

The kernel receives contributions from Sa ,/3 and r^, and can be written as 
i 7 |(fci, k2\a) = -7(a) [G^(fci, ^2) + Ga^{ki, ^2)] - 7 i(«)G'i(fci, ^2) - 72(0)02(^1, ^2) 
where we have defined 


7 i(a) = - 


1 


da' Gs{a, a') a'gi{a') 


(4.38) 

(4.39) 

(4.40) 


Explicit expressions for the kernels G^^^ai3,i} a-i'e given in Appendix G. Let us note that the value 
of 7 ( 0 ) has already been fixed by the renormalization of the power spectrum and therefore can no 
longer be adjusted. Substituting (4.34) for 7 (a), the sum of the loop diagram and the counterterm 
becomes 


-^113 A Bile 








— “/n 


/ 74 
I - r 


f 52 
I - f 


^^^g^cT 2 (A) + 71 (a) j kl + ^^a^(A) + 72 ( 0 ) 
X Pi-^{ki)Pii{k 2 ) + 5 perms . 


(fc3 • k2 

kl 


- -ki 


(4.41) 
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The cutoff dependence is therefore cancelled if 


71 (A, a) 

72 (A, a) 


74 

-^cr2(A, a) + finite 77 
5o7 

52 

—cr^(A, a) + finite 77 
loo 


51 (A, a) 

52 (A, a) 


74 1 2/. N . . 

(A,o) Thnite , 

DO a 
52 1 

— —-(T^(A, a) + hnite . 


(4.42) 

(4.43) 


We have shown that every counterterm generated by the viscosity contribution to the stress 
tensor, r^, yields the correct momentum dependence to remove some of the loop divergences 
in the power spectrum and the bispectrum. Next, we show that the remaining divergences are 
cancelled by the noise counterterms. This part is for aficionados and can be skipped on a hrst 
reading. 


Noise counterterms 

Noise terms appear both in the initial conditions (3.8) and in the stress tensor (3.35). They are 
essential in the renormalization procedure as they remove some divergences in P 12 and as 

well as all the divergences in B ^22 • 

As indicated in (4.25), the noise term can be expanded in powers of the long-wavelength 
fluctuations. The n-th order solution <5/ n receives contributions from both Gaussian and non- 

(n) 

Gaussian initial conditions: 


Specihcally, at one-loop order, we will require the first-order Gaussian contribution 

6'^^[^{k,a) = kikjNl\k,a) , (4.45) 

and the first and second-order non-Gaussian contributions 
6l[^^{k,a) = UM{k)',Pj{k) , (4.46) 

A7s^^(A:,a) =/nl / hkjNlij{k - p,a)ij{p) + 2F2{k - p,p)6(^i){k - p,a)M{p)'tljj{p) , (4.47) 

where 


s-J,NG/ 


(4.44) 


Gs{a,a)jQ^^^{k,a ) . 


(4.48) 


We will now show that these connterterms cancel the divergences proportional to M(k) in both 
P 12 and and all the divergences in 3^22 ■ 

• The divergence proportional to M{k) in P 12 is renormalized by the counterterm More 
precisely, we have 

P 12 + PJJ 3 -l;Ual{K)k^M{k) + 2{5l{^^{k,a)5^,\^{-k,a))' 


21 

= /nl 


hi) 


-^cr^(A) + 2 / da' (a, o') Jo (a) 

J din 


k^M{k) = finite , (4.49) 
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where Jo{a') was defined in (3.38). Hence, the cutoff dependence is eliminated provided 
that 

^o(A, a) = a) + finite . (4.50) 

ba ^ 

Similarly, the divergence proportional to M(k) in B ^22 (which is the same as the divergence 
in P 12 ), is renormalized by the second term in (4.47). 

• Let us consider the following correlation between the first-order SPT solution and the noise 
terms 


Buj 3 + (5fPT(fci, a)h)'()^(fc 2 , a))' 

^ B[fj + B[fj (4.51) 

The first term, B^^j, cancels the divergence proportional to M{k) in while the second 
term, h|jj, cancels the first divergence in B^ 22 - demonstrate the second fact 

explicitly and leave the first as an exercise to the reader. 

Using (4.46) and (4.47), we can write 

BiTj = {k 2 )iik 2 )j{k 3 )k{k 3 )i{N"J{k 2 ,a)N!^''{-k 2 ,a)y Pi^{ki) + 5 perms , (4.52) 


where was defined in (4.48). We will assume that the cutoff-dependent part of 
satishes oc a, so that — and 


(iV*^(fc,aX(-fc,a))' = -. 
Hence, we find that 


J{a)5^^5^^ + J{a) ( 5 ^^ 5 ^^ + S 


(4.53) 


?(i) 




-°122 w Rijj Jnl 






Pi'ipiki) + perms 


= finite . (4.54) 

We see that b[^j has the right momentum dependence to absorb the hrst divergence in H ^22 • 
In particular, the cutoff dependence is cancelled if 

13 - 20 1 

J'(A,a) = —^d'2(A, a)-|-finite and J(A,a) = -^it 2(A, a)-|-finite . (4.55) 

3a^ 3 

• Finally, we consider 

Bjjj D a)(5)^’|^(^ 2 , (^ 3 , a))' + perms 

= /NLM(A:i)(A:2)*(fc2),(A:3)fc(fc3)z(V^j(fci)JS''(fc2,a)4'(^3,a))' . (4.56) 


24 






It is difficult to determine the precise momentum dependence of this correlation function, 
since the probability distributions of the noise terms are unknown. We will assume that it 
can be expanded around fc* = 0, so that we can write 

{i:j{ki)Jl\k2,a)4\k^,a))'=M{a)5^^8'^^+N{a){5^'^6^^ + 5^5^'^)+0{kl) . (4.57) 


In that case, we get 


-^122 + ^ /n 


+ 


+AA(a) ) klkl 


80 

147 


+ 2AA(a) ) (fc2 • ks) 


M{ki) + perms 


= finite . 


(4.58) 


We see that Bjjj has the correct momentum dependence to absorb the cutoff dependence 
of the second divergence in (4.15). More precisely, we have 
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AA(A, a) = —(T^(A, a) + finite 


and M{A,a) = — y^< 5'^(A, a) + finite . (4.59) 


We have shown that the solutions generated by the stress tensor derived in Section 3 have the 
correct momentum dependence to cancel all divergences coming from SPT loop diagrams. So 
far, we have focused on the cutoff-dependent parts of the EFT parameters. In the next section, 
we study their finite (or renormalized) parts. 

4.3 Renormalized EFT Parameters 

Any EFT parameter can be written as the sum of a cutoff-dependent part (the counterterm) and 
a finite (or “renormalized”) part: 

g{A, a) = 5 ^°) (A, a) -h (a) . (4.60) 

Ultimately, every long-wavelength observable must be independent of the cutoff. Indeed, once 
the large-scale correlation functions have been properly renormalized, one can send the cutoff to 
inhnity. In ACDM, all loop integrals are convergent and the cutoff can be taken to infinity even 
before renormalization. This may seem to go against the EFT philosophy since we are including 
modes which are outside the regime of validity of the EFT. However, the mistakes one makes in 
doing so, can be absorbed into a shift of the renormalized EFT parameters. 

The time dependence of the renormalized parameters may be different from that of the coun¬ 
terterms. Moreover, it cannot be computed within the EFT framework, but must be determined 
from simulations or observations. An exception is a “scaling universe” — i.e. a matter-dominated 
Einstein-de Sitter (EdS) cosmology with scale-free initial conditions Pn{k, ain) oc kA. In this spe¬ 
cial case, an additional symmetry constrains the time dependence of the parameters in the EFT. 
To see this, we first note that the equations of motion in EdS are invariant under a Lifshitz sealing 

X I— )• XxX and a i—)• XaO . (4-61) 
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The dimensionless power spectrum of the initial conditions, oc is also invari¬ 

ant under the Lifshitz scaling iff 

Aa; = . (4.62) 

When (4.62) holds, then the Lifshitz scaling maps one EdS solution to another with a different 
realization of the same statistical initial condition. The linearly-evolved power spectrum can 
be written in a manifestly self-similar form, A^]^(A:,a) = (A://cnl)"^"*~^) where the nonlinear scale 
satisfies /cnl oc Nonlinear corrections to Af^ appear as higher powers of /c/Zcnl [6, 37]. 

Assuming that the short-scale fluctuations satisfy the Lifshitz scaling, the terms in the effective 
stress tensor must also have the right transformation properties. In particular, the symmetry 
fixes the combination of powers of k and a which can appear in the equations of motion and 
therefore constrains the time dependence of the renormalized parameters. For instance, using 
such an argument,® one finds that the speed of sound in a scaling universe must satisfy c^(a) oc 

For non-Gaussian initial conditions, we also need to impose that the higher-point correlation 
functions are invariant under the Lifshitz scaling (4.61) with A^, and Aa related by (4.62). More 
specifically, let us assume that the dimensionless bispectrum satisfies 


Bin 


f k^ k^ k^ 

\ Ax Ax Ax 



X^XlBiiiikiMM^ain) 


(4.63) 


where we have used that Bni{ki-,ain) oc Given (4.62), the initial bispectrum is invariant 
under (4.61) iff m = — |(n -|- 3). Furthermore, if the primordial Gaussian potential (pg is scale- 
invariant, then the field ip must transform as '!/’(*) '4’i^xx) under the Lifshitz scaling; 

cf. eq. (3.13). Using that the Lagrangian coordinate transforms as q{x,a) i—)• Xp^q{XxX, Xao), we 
find that the field ^{x,a) = 'ip{q{x,a)) transforms as 


T(a;,a) A^T(Ax®,Aaa) . 


(4.64) 


This implies that the terms in the stress tensor (3.31) preserve the Lifshitz scaling iff their time 
dependence is 

, gipia)} oc oc . (4.65) 

In AGDM and for general initial conditions, the Lifshitz scaling isn’t a symmetry and the initial 
conditions aren’t scale-free. However, in fig. 4 we show that the initial matter power spectrum has 
momentum regimes in which a power law ansatz is a good approximation. This explains why, for 
Gaussian initial conditions, the ansatz deduced from the Lifshitz scaling symmetry has worked 
remarkably well [8, 9, 36]. We will extend these results and assume that the EFT parameters 
have the following time dependence^ 

g(R){a) = [Hia)fia)f[Di{a)ro+^g^^^ , (4.66) 

®Let us remark that this argument does not determine the time dependence of the counterterms a), since 

the cutoff introduces a new scale and therefore breaks the Lifshitz scaling symmetry [6]. 

^Notice that we have included a factor {Hf)^ which is not present in [8]. We found that this ansatz improves 
the accuracy of (C.IO) by a few percents. 
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fc[/iMpc 

Figure 4: Scaling of the initial matter power spectrum. We see that in the momentum regions 
[0.02, 0.07] /iMpc“^ and [0.07,0.25] /iMpc“^ the power spectrum is well approximated by a power 
law with n « —0.9 and /cnl ~ 0.16/iMpc“^ (red line) and n « —1.5 and ~ 0.23/iMpc“^ 
(blue line), respectively. 


where mg is a real parameter and is a constant. At one loop, we don’t need to know the 
precise time dependence of gl 2 to compute the counterterms. The ansatz (4.66) is therefore only 
required for the parameter g{a). In general, the value of mg is unknown. In what follows we 
will assume that the time dependence matches the ansatz of the scaling universe, i.e. we will use 
rUg = (1—n+2A)/(n+3). From fig. 4 we deduce that n G [—1.5,—0.9] for/c G [0.02,0.25]/iMpc“^. 
For these values of n, we have mg G [0.9 + A, 1.7 + 1.3A]. For our numerical computations, we 
will choose mg = 1.7 + 1.3A. However, we have checked that the counterterms only change by 
at most a few percents as the value of mg runs over this interval, so our numerical results will 
not depend sensitively on this choice. 

5 Numerical Analysis 

In this section, we present a numerical computation of the renormalized power spectrum and 
bispectrum for non-Gaussian initial conditions. We start, in §5.1, with a summary of the different 
contributions (SPT and EFT) to the dark matter correlation functions. In §5.2, we estimate the 
relative sizes of the terms in a scaling universe. Finally, in §5.3, we make some preliminary 
observations on the shapes of the various contributions. A more detailed analysis will appear 
in [18]. 

5.1 Renormalized Correlation Functions 

In §3.3, we wrote the effective stress tensor of the dark matter as a sum of terms made from the 
long-wavelength fields and their derivatives 

Tv — ^ ^ Qn Qn ) (^'l) 

n 
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where the coefficients an and the corresponding operators On relevant for this paper are listed in 
Table 1. 


On 

A 5 

A(52) 

A(s2) 

d\sijd^ 5 ) 

AT - di{ 5 d^^) 

A(T(5) 


dn 


ei 

62 

63 

9 

91 

92 

CXn 


61 

62 

63 

7 

71 

72 


Table 1: Operators and parameters in the EFT-of-LSS with primordial non-Gaussianities. 


In §4.2, we showed that these operators gives rise to counterterms which renormalize the one-loop 
power spectrum and bispectrum of standard perturbation theory. Schematically, the solution can 
be written as (see §4.2 and §C.2 for more details) 


S^{a) ~ 



da Gs{a,a)Tn{a') ~ '^an{a)Onia) , 

n 


(5.2) 


where the parameters an can be expressed as integrals over time of the corresponding parame¬ 
ters an- The SPT solution (a) and the EFT counterterms 5‘^{a) generate correlation func¬ 
tions for the renormalized dark matter density contrast 6 (a) which depend both on the cos¬ 
mological parameters of the ACDM model 9 = /i, A<^, n^} and on the EFT parame¬ 

ters a = {^, 61 , 2 , 3 , 7 , 71 . 2 }- 


• The dark matter power spectrum can be written as P = + /nl-P^*^, where P^ arises 

from the nonlinear evolution of Gaussian initial conditions [5] and P^^ is the contribution 
from non-Gaussian initial conditions. Both P^ and P^^ can be split into an SPT part and 
an EFT part. At one loop, we have 

P^T = -^11 + -^13 + P 22 , Peft = —, (5.3) 

Pg'^PT = P 12 , Pef^'t = , (5.4) 


where Pi^(/c) was defined in (4.2). 


• The total dark matter bispectrum is P = + /nlP^*^, where is the part arising from 

the nonlinear evolution of Gaussian initial conditions [8, 9], and B^^ is the part associated 
with non-Gaussian initial conditions. Both B^ and B^^ can be split into an SPT part and 
an EFT part. At one loop, we have 


pG _ 

-^SPT “ 


P 112 + 


P 114 + P 123 + P 123 p P 222 


(5.5) 


PsPT ~ Pm + 


B 


(I) 

113 


d(II) I nw I TJ 
-Dill w -£->109 w 


Tl3 


(I) 


122 


9 ( 11 ) 

T 22 


(5.6) 


3 

Peft ~ P , (5-7) 

i=l 

2 

Peft = ^Pf P ^^7 P E P^7, • (5-8) 

i=l 
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The individual bispectra in (5.7) and (5.8) can be written as 


Bf =-2[Eai3iki,k2) + Es{ki,k 2 )]Pii{ki)Pii{k 2 )+ 2 perms , (5.9) 

^NG ^ kl)Bni{ki,k2, fcs) , (5.10) 

B^^ = -2Ei{ki,k2)Pii{ki)Pii{k2) + 2 perms , (5.11) 

B^ = - [Gai3iki,k2) + Gii,{ki,k 2 )]Pii{ki)Pi^{k 2 ) + 5 perms , (5.12) 

B^^ = -Gi{ki,k 2 )Pii{ki)Pi^{k 2 ) + 5 perms , (5.13) 


where the kernel functions E..{ki, k 2 ) and G..{ki, k 2 ) are defined explicitly in Appendix C. 

For simplicity, we will take the cosmological parameters to be fixed and allow only the EFT 
parameters to vary (for a more complete treatment see [18]). Furthermore, we will assume 
that the sound speed parameter has been measured in the power spectrum and its value is 
fixed ^ = 1.5 h~‘^ Mpc^ [ 8 ]. The following parts of both the Gaussian and the non-Gaussian 
bispectra are therefore predicted 



Bl 

= Bspt + ^B^ 

I = G,NG. 

(5.14) 

The total bispectrum can 

then 

be written as 




B = 

B^ + B^ + /nl 

/rNG , PNG'! 

[Bq -I- j 

, (5.15) 

where we have defined the 

sum 

of the additional counterterm contributions as 



i=l 


(5.16) 



S^G = + 

2 

(5.17) 


i=l 


In general, the coefficients in (5.16) and (5.17) cannot be predicted, but need to be measured 
in N-body simulations or in observations. For the numerical results of §5.3, we will estimate 
the sizes of these coefficients by looking at the one-loop divergences that they cancel. In 
particular, since in our universe loop integrals are convergent, they are in practice evaluated 
with the cutoff A taken to infinity. In that case, the SPT loop diagrams contain integrals 
of the linear power spectrum extrapolated to scales beyond the linear regime (i.e. k > /snl)- 
This is not a problem, since, as we explained in §4.3, this finite error is removed by adjusting 
the counterterms. More precisely, from (4.34) we infer that the term in 7 which cancels 
this known error is 


I7I ~ 


2 



dp gp 
2vr2 (p/p)^ 


Pu{p)-\ - 


(5.18) 


where /cnl is the nonlinear scale. We will choose /cnl = 0.2/iMpc“^, which corresponds to 
the nonlinear scale found in [ 8 , 9]. For local non-Gaussianity (i.e. A = 0), we then find 
I 7 I ~ 1.6 h~^Mpc^. Note that, for A 0, the precise value of p is unimportant as it cancels 
in the product 'yB^. 
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Of course, there is also an unknown, finite contribution to 7 that accounts for the effects 
of the nonperturbative short scales on the long-wavelength dynamics. We will assume that 
there is no fine-tuned cancellation between this finite part and the term shown in (5.18). 
The result in (5.18) then gives an approximate lower bound on the size of the counterterm. 
Applying the same estimate to 71 and 72 , we find 

71 ~ ^7 and 72 —7 . (5.19) 

Let us stress that the estimates in (5.18) and (5.19) only provide approximated lower bounds 
and should not be considered precise evaluations. 


In §5.3, we will present numerical results for three different primordial bispectrum shapes: 


Local non-Gaussianity is obtained by setting = 1 in (2.18). The primordial bispectrum 
then is 

■ 1 1 1 


+ 


+ 


lSlS 1.3^3 1-3 j,3 ’ 

/vi a^2 ^1*^3 ^2^3 


(5.20) 


where A = 27r^A^. For simplicity, we have written the bispectrum for scale-invariant initial 
conditions, ns = 1. The latest CMB constraint on the amplitude of the local bispectrum is 
= 1.8 ± 5.6 [38]. 


• Higher-derivative corrections to slow-roll inflation produce equilateral non-Gaussianity [24, 
25]. In the effective field theory of inflation [25], this is captured by two cubic operators, 
TT^ and 7r(9j7r)^, for the Goldstone boson of broken time translations, vr. Both operators 
produce very similar equilateral bispectra. For purposes of illustration, we focus on the 
PNG produced by the operator vr^, which has a bispectrum of the form 

fe) = 162 /Sr ■ 5)5 . (5.21) 

where K = ki + k 2 + k^. The latest GMB constraint on the amplitude of the equilateral 
bispectrum is = —9.2 ± 69 [38]. 


• The shape of the bispectrum in quasi-single-field inflation cannot be computed analytically. 
However, it is well approximated by the following ansatz [26, 39] 




1 Nu[8k] 

7^N48/27] ’ 


(5.22) 


where k = kik 2 ks/K^ and Ni, is the Neumann function of order u. The parameter u 
depends on the mass of the hidden sector held to which the inhaton held couples during 
inhation. It determines the scaling of the bispectrum in the squeezed limit via A = ^ — u 
in (2.20). For purposes of illustration, we choose 5 corresponding to A = 1 which is 
precisely intermediate between the scaling of local and equilateral non-Gaussianity. 


In our numerical computations, we consider a hat ACDM cosmology with the following pa¬ 
rameters: = 0.27, = 0.73, h = 0.70. The amplitude of the primordial potential is 
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A^(A:o) — 8.7 x 10“^*^, defined at the pivot scale = 0.002/iMpc“^. The one-loop correlators 
are computed numerically with Mathematica. For validation purposes, we wrote two independent 
codes. Both codes yielded results which are in very good agreement.® 

Our results are summarized in figs. 5, 6 and 7 for local, equilateral and quasi-single-field 
PNG, respectively. The plots show the SPT and counterterm contributions to the bispectrum 
separated into Gaussian terms (which are the same in every plot) and non-Gaussian terms. 
Here, and in the following, the results are evaluated at redshift z = 0. Notice that while PNG 
has a larger effect on larger scales, we have chosen to show the “mildly nonlinear” range of 
scales, where the EFT counterterms can play an important role. For each type of PNG, we 
represent the corresponding matter bispectrum in the equilateral configuration (left) and for 
fixed ki = 0.01/iMpc“^ (right). We see that, in the equilateral configuration, the two-loop 
Gaussian contribution becomes comparable to the non-Gaussian contribution on rather large 
scales. In §5.2, we provide an analytical estimate of the critical scale kc at which the two-loop 
Gaussian contribution becomes relevant. On the other hand, we see that the two-loop contribution 
is much more subdominant in squeezed configurations. This is also the momentum configuration 
in which the shape of PNG (and more precisely, the scaling A of the squeezed limit) can leave an 
appreciable imprint on the matter bispectrum. In §5.3, we discuss the possibility of extracting 
information about the shape of PNG from the matter bispectrum. 




fc[ftMpc 


Figure 5: Gontributions to the bispectrum for local non-Gaussianity with = 10, evaluated 
in the equilateral configuration {left) and for fixed kL = O.Ol/iMpc”^ {right). The definitions of 
and can be found in (5.14), (5.16) and (5.17). 


5.2 Estimates in a Scaling Universe 

In this section, we determine the critical scale kc at which the two-loop Gaussian contribution 
becomes relevant. For purpose of illustration, we will consider local PNG, although our result 

® As a technical aside, let us note that the expression of the one-loop bispectrum in (5.6) is ill-suited for numerical 
evaluation as the different contributions each contain IR divergences. The equivalence principle guarantees that 
these divergences cancel when all diagrams are summed [37]. However, to avoid cancelling large numbers against 
each other, it is preferable to rewrite the integrand in a form which makes this cancellation manifest. In Appendix D, 
we present the “IR-safe” integrands for the non-Gaussian loops. 
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Figure 6: Same as fig. 5 but for equilateral non-Gaussianity with = 100. 



won’t depend sensitively on the precise shape of the non-Gaussianity. To estimate the relative 
sizes of the various contributions to the power spectrum and bispectrum we consider a scaling 
universe (i.e. EdS with power-law spectrum Ps{k,Tin) oc fc"'). This extends the results of [6, 8] to 
non-Gaussian initial conditions. 

As explained in §4.3, in a scaling nniverse, Lifshitz transformations of the form (4.62) leave the 
initial dimensionless power spectrnm unchanged, which can therefore be written as a polynomial 
in /c/Zcnl- Because the evolution is self-similar, the same applies to the nonlinear power spectrnm 

= ^^sik,T) = A|(/c/A:nl) • (5.23) 

Furthermore, the transformation (4.63) implies that the initial bispectrum is also self-similar, 
and that the nonlinear dimensionless bispectrum is therefore a polynomial in ki/k^i,: 

f k^ 

Bsiki,k2,k3,T) = Bs{ki,k2,k3,T) = Bs{ki/kj^j^,k2/kt^i^,k3/kj^i^) . (5.24) 

The different contributions to the dimensionless power spectrum and bispectrum will scale as 
different powers of k/k^i^. 
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Power Spectrum 

Let us first collect the terms that contribute to the power spectrum A|(A;/fcNL)- The momentum 
scalings of the Gaussian contributions were derived in [6]: 




C 


s{n) 



p 


n -|- 3 ? 

a 1.5 

tree 

2{n + 3) f 

3.0 

1-loop 

n -|- 5 ? 

3.5 

LO vis. 

3(n-F3) ? 

a 4.5 

2-loop 

2n -|- 8 f 

5.0 

NLO vis. 

n -|- 7 f 

5.5 

NLO h.d. 

4(n -I- 3) ? 

a 6.0 

3-loop 

7 f 

7.0 

noise 


where ‘tree’ refers to the linear power spectrum ‘loop’ stands for the SPT loop contributions, 
‘LO vis.’ is the leading-order viscosity counterterm defined in (5.3), ‘NLO vis.’ are the next-to- 
leading-order viscosity counterterms [6], ‘NLO h.d.’ are next-to-leading-order higher-derivative 
terms, and ‘noise’ refers to Ajj. The overall coefficient s(n) is a (computable) order-one number, 
which is different for the different contributions (see Appendix A of [6]). For the numerical 
estimates of the scaling index p we have used n ~ —1.5, which corresponds to the scaling of 
the linear power spectrum in the regime k G [0.07,0.25]/iMpc“^ (see fig. 4). The scaling for 
k G [0.02,0.07] ZiMpc”^, which is n ~ —0.9, will also be relevant, but hasn’t been shown explicitly. 

For the non-Gaussian contributions to the power spectrum, we find 




c 


s{n) /nlA,^ 



P = 

= i(^ 

+ 3) 

a 2.25 

1-loop 

P = 

= 

+ 7) + A ? 

a 2.75 + A 

LO vis. 

P = 

= i(^ 

+ 3) 

Si 3.75 

2-loop 

P = 

= 

+ 3) 

a 5.25 

3-loop 

P = 

= 

+ 13) ? 

a 5.75 

noise 


where ‘LO vis.’ refers to the leading order non-Gaussian viscosity counterterm, defined in (5.4), 
‘loop’ stands for the SPT loop contributions and ‘noise’ stands for Ajj. The scalings of the 
viscosity counterterm and the noise term are fixed by the /^-dependent parts of the UV-limit of 
Af 2 (see §4.1). We note that the amplitude of the non-Gaussian contributions is suppressed by 
a factor of /nlA<^ ~ 3 x 10“^/nl- Higher-loop corrections are suppressed by additional factors of 
Aii(^) = (^/^nl)*^"*"^- As before, the overall coefficient s{n) is a computable order-one number 
which is different for the different contributions. 
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We see that the non-Gaussian contributions to the power spectrum are highly suppressed. 
We should therefore ask when the Gaussian two-loop terms (which we haven’t included in our 
analysis) are of the same order. The two-loop Gaussian contribution is estimated through a single 
representative term, namely Pg® [7]: 

^ 2 Uoop C pS = ^ ^ (5.25) 


A rough estimate of the scale kc at which the two-loop Gaussian contribution equals the one-loop 
non-Gaussian contribution is 


pG 

^ 2-loop 

pNG 

'^l-loop 


/nlA(^ Si2 V/Cnl/ 


0.18 A)nl 



(5.26) 


where S 33 and S 12 are the coefficients of P 3 ® and P 12 , respectively. For the numerical estimate 
in (5.26) we have used® n k, —0.9 , /cnl ~ 0.16/iMpc“^ and S33/S12 — 0.06. In hg. 8, we compare 
this estimate for kc{/ nl) with the value obtained from a numerical computation in a ACDM 
universe. We see that already on relatively large scales the two-loop Gaussian term is of the same 
order as the one-loop non-Gaussian term. 




Figure 8: Comparison between the 1-loop contribution to the power spectrum for local PNG 
with /nl £ [1,100] (grey band) and an estimate of the 2-loop Gaussian contribution (red dashed 
line). The right plot shows the critical momentum value kc at which the two-loop Gaussian 
contribution becomes larger than the one-loop non-Gaussian contribution as a function of /nl, 
evaluated numerically (solid line) and using the estimate (5.26) (dashed line). The numerical 
result is for a ACDM cosmology, while the estimate is for a scaling universe. The shaded region 
represents the momentum range where the scaling n ~ —0.9 is no longer valid in the real universe. 


®These values are applicable to the momentum range k £ [0.02,0.07]/iMpc“^ (see fig. 4), which is roughly 
where the two-loop Gaussian contribution is expected to be larger than the one-loop non-Gaussian contribution 
(see fig. 8). 
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Bispectrum 


For the bispectrum, we focus our attention on the equilateral configuration, Bs{k, k, k), for which 
the two-loop corrections are largest. For Gaussian initial conditions, the scalings of the various 
contributions were derived in [ 8 ]: 


Bf{k) C s(n) 


2(n-h3) ? 

3.0 

tree 

3{n + 3) ? 

=i4.5 

1 -loop 

2{n + 4) ? 
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NNLO vis. 

2 n-hl 0 ? 

a 7.0 

NNLO h.d. 

5(n-h3) ? 

7.5 

3-loop 

n -|- 10 ? 

a 8.5 

noise 


where ‘NLO vis.’ and ‘NNLO vis.’ refer to the next-to-leading-order and next-to-next-to-leading- 
order contributions to the Gaussian counterterm, respectively. More precisely, the ‘NLO vis.’ 
term corresponds to the term defined in (5.7), while the ‘NNLO vis.’ term corresponds to 

the one-loop diagrams formed with the leading Gaussian counterterms. The term ‘NNLO h.d.’ 
represents higher-derivative contributions and the term ‘noise’ corresponds to Buj. 

For the non-Gaussian contributions to the bispectrum, we find 


B^ik) 
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s{n) /nlA^ 
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where ‘LO vis.’ and ‘NLO vis.’ refer to the leading-order and next-to-leading-order contributions 
to the non-Gaussian counterterm respectively. More precisely, the term ‘LO vis. (1)’ refers 
to B^® (5.10), the term ‘LO vis. (2)’ refers to B^*^ (5.17), and the term ‘NLO vis.’ refers to the 
finite part of the one-loop diagrams formed with one counterterm. The term ‘NLO h.d.’ represents 
higher-derivative contributions and the terms ‘noise (A)’, ‘noise (B)’ and ‘noise (G)’ refer to the 
noise terms b[jj, B^jj and Bjjj (cf. §4.2). As before, the scalings of the viscosity and noise 
contributions to the bispectrum are determined by the fe-dependence of the UV-limit of the loop 
diagrams they cancel (cf. §4.1). 
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Figure 9: Top: Comparison between the tree-level contribution to the bispectrum for local PNG 
with /nl £ [1) 100] (grey band) and an estimate of the 2-loop Gaussian contribution (red dashed 
line). Bottom: Comparison between the 1-loop non-Gaussian contribution (gray band) and the 
2-loop Gaussian contribution (red dashed line). The curves in the right plots have the same 
meaning as in fig. 8. 


Since the non-Gaussian contributions to the bispectrum already appear at tree level, they are 
less suppressed than in the power spectrum. Again, we should ask up to which point we can 
neglect the Gaussian two-loop corrections. To estimate the two-loop corrections, we will consider 
one representative contribution, namely PI- 

-®&oop 'T .0233 ~ 

We can then determine the scales at which the two-loop bispectrum starts to become important 


(5.27) 
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relative to tree-level bispectrum Bm and the non-Gaussian one-loop contribution 

(5.28) 

(5.29) 

As before, we have used n ~ —0.9 to estimate the numerical values of kcif^i.)- For the numerical 
prefactors we have used — 0-11 and •S 233 /sf_)^op — 0.01. In fig. 9 these estimates are 

compared to the exact result of a numerical computation in ACDM. Again, we find that two-loop 
corrections becomes relevant at relatively large scales. 

We wish to emphasize that the importance of Gaussian two-loop corrections is maximal in the 
equilateral configuration. Away from this limit the two-loop amplitude is significantly smaller. 
In hg. 10, we show tree-level and one-loop contributions associated with local non-Gaussianity 
for fixed kL = 0.01 /iMpc“^. This time the two-loop correction is subdominant for all scales, even 
for relatively small values of /nl- 



Figure 10: Comparison of the non-Gaussian tree-level and one-loop contributions with an esti¬ 
mate of the Gaussian two-loop contribution, evaluated for fixed k^ = 0.01 fiMpc”^. 
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5.3 Comments on Primordial Shapes 

The squeezed limit of the bispectrum contains interesting information about the spectrum of 
particles during inflation. In particular, it can reveal the masses [26, 40] and spins [29] of particles 
to which the inflaton couples. This requires an accurate measurement of the scaling dimension A 
in the primordial bispectrum (2.20). In [39], the detectability of this signature was discussed for 
future galaxy surveys. However, their analysis only used the tree-level form of the dark matter 
bispectrum, so we should ask when the nonlinear corrections discussed in this paper become 
important. (See [41] for related observations in the special case of local non-Gaussianity.) 
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Figure 11: Comparison of the non-Gaussian SPT contribution Bgprj, for local, equilateral and 
QSF shape, normalized to the local tree-level bispectrum Dashed lines refer to non- 

Gaussian tree-level contributions, while solid lines include one-loop corrections. We see that the 
loop corrections become relevant on relatively large scales and tend to decrease the difference 
among the three primordial signals. 

In fig. 11, we show the non-Gaussian contributions to the bispectrum normalized with 

respect to the tree-level bispectrum of the local type We have plotted three different 

primordial bispectra — local, equilateral and quasi-single-field — whose shapes were defined 
in §5.1. The dashed lines show the different scaling behavior of the tree-level contributions. We 
see clearly that the one-loop corrections become important at relatively large scales, especially as 
we go further away from the scaling of local PNG. Moreover, the loop contributions are much less 
sensitive to the value of A. In fact, the one-loop contributions for equilateral and quasi-single- 
field scale almost the same. This suggests that adding the one-loop corrections is important and 
that it can make it harder to extract A (relative to the tree-level expectation). We will study 
this question in more detail in [18]. 

6 Conclusions 

Measurements of primordial non-Gaussianity in the cosmic microwave background are nearly 
saturated — correlations on large angular scales can’t be measured more precisely because of 
cosmic variance, while Silk damping limits the information that can be extracted from small 
angular scales. In contrast, constraints on PNG from LSS surveys are only starting to become 
available and have the potential to improve significantly in the future [1]. However, to extract this 
information requires pushing the theoretical understanding of gravitational clustering to smaller 
scales into the mildly nonlinear regime. 

Nonlinear gravitational evolution produces two types of effects: First, it generates non- 
Gaussianity even if the initial conditions are perfectly Gaussian. This inevitable “background” 
needs to be characterized precisely before any non-Gaussian “signal” can be discovered. Second, 
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any primordial non-Gaussianity in the initial conditions gets distorted. These distortions can 
affect how well we can distinguish between distinct types of non-Gaussianity (cf. fig. 11). Both of 
these effects are captured by the EFT-of-LSS. The first was studied in [8, 9], while, in this work, 
we have developed the tools to address the second. 

A number of open questions have been left for future work [18]: 

• For small levels of PNG, two-loop Gaussian corrections become comparable to the one-loop 
non-Gaussian corrections already on relatively large scales. (This is especially true for the 
power spectrum and the equilateral limit of the bispectrum.) Extending the theoretical 
treatment to fc ~ O.l/iMpc”^ (and beyond), therefore requires a more complete treatment 
of the two-loop corrections for Gaussian initial conditions. 

• The hnite parts of the EFT contributions cannot be predicted, but have to be measured 
in N-body simulations or in observations. For Gaussian initial conditions this one-loop 
matching has been performed in [5, 8, 9]. The one-loop matching for non-Gaussian initial 
conditions is still outstanding. 

• Naively, the number of useful modes in LSS scales as the cube of the maximum wavenumber, 
kmax, at which theory errors are still under control. However, in practice, nonlinear evolu¬ 
tion moves information from low-order correlations to higher-order correlations [42, 43]. It 
remains to be quantified how much information can actually be extracted from measure¬ 
ments of the dark matter bispectrum. 

• Most LSS observations only provide access to the statistics of biased tracers (galaxies, halos, 
etc.) of the underlying dark matter density field. The biasing between these tracers and 
the dark matter introduces an additional source of nonlinearity. On the other hand, the 
scale-dependence of the biasing also provides the opportunity of seeing specihc imprints of 
the non-Gaussian initial conditions [44]. A consistent model of nonlinear biasing requires 
renormalization [10-13]. The first steps towards defining a self-consistent nonlinear biasing 
model in the presence of non-Gaussian initial conditions have appeared in [15]. 
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A Odd Spin in the Squeezed Bispectrum 


In this appendix, we make a few comments about the allowed angular dependence in the squeezed 
limit of the bispectrum (see also [28, 45]). We will explain why, at leading order, we should expect 
only even-spin contributions to the squeezed limit. 

The following two types of squeezed limits are commonly used in the literature: 


Wm. B^p{\\k + p\,\\k - p\,k) , (A.l) 

lim B^{p', \p' -k\,k) . (A.2) 

fc—>-0 


In this work, we have defined the Legendre expansion of (A.l), 


lim B^{\\k + p\,\\k - p\,k) = 




. L.i 


p 


PUk-p) 


1 + 01 ^ 

pz 


P^{p)P^{k) . (A.3) 


Since the left-hand side is invariant under p i—?■ —p (regardless of whether p is large or not), we 
must have ai^i = 0 for L odd. Comparing (A.l) and (A.2), we find that p' = p + and hence 


k-p = k-p' - -[l -{k- p'Yj + O 


{p'f 


(A.4) 


Substituting (A.4) into (A.3), we get 


k 


Phik .p) = PL{k- p') + ^(-^ + 1) \ PL+i{k • p') - {k • p')PL{k • p') 


(A.5) 


This shows that we do generate Legendre polynomials of odd powers in the squeezed limit (A.2). 
However, was also see that they are always subleading since they are forced to come together 
with a leading Legendre polynomial of even degree. 
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B Coarse Graining in Perturbation Theory 


In this appendix, we explicitly integrate out perturbative short scales in the Euler equation and 
show that all terms consistent with the symmetries are generated in the resulting effective theory. 

Suppose that the fluid equations are defined at a scale A < and we wish to integrate out 
modes in the momentum shell p G [A, A], The result is an additional stress tensor in the Euler 
equation [4]: 

+ [pvlvi]^ . (B.l) 

We wish to expand in terms of the long-wavelength fields. We will illustrate the computation 
for a single term 

. (B.2) 

For simplicity, we will work in an Einstein-de Sitter universe and compute to first order in 
the long-wavelength fluctuations. 


B.l Nonlinear Evolution 


Since all modes are perturbative, can be computed in perturbation theory. Below we show 
that 


a’‘^{k,T) = 


d^p 


JA (27r)3 
where we have defined 


rA 


p^p^ap{k) + 


d^k 

W) 


Si^{k,k)ap{k-k)6{k,T) + --- 


3 P 


(B.3) 


ap{k) = [(l)s{p+ \k,Tin)(t)s{-P + \k,Tin)\f^ , 


s;^ik,k)=pY 


20 ^ {k — k) ■ k 


20 {p ■ k)'^ 

7 p2'f,2 




+ 

P*p +p^k 


p ■ k 

~w 


(B.4) 

(B.5) 


The first term in (B.3) simply arises from the linear evolution. The second term comes from 
replacing one of the (p by its second-order solution (/)( 2 ) and extracting the dependence on the 
long-wavelength fluctuations. We give the details of the derivation in the following insert. 


Derivation. The first-order contribution to (B.2) is 


— = 
97^41 Jfi) 


( 1 ) 


(^ 


{p,k) 5i{-p + \k)52{p + \k) + 52{-p + \k)5i{p + \k) , (B.6) 


where we have used —kPpik) = ^'H^5(k) and defined 


ff'(p,fc) = 


+ W) 

{p- \kY{p+ lfc)2 


(B.7) 


Substituting (C.15) for ^ 2 , we get 

d^p 


9774 


= / n*4(p,fc)7^2(fe,p+5fc-fc)5i(-p-t-|fe)i5i(p-b ifc-fc)(5i(fe) . (B.8) 

^ ' JA Jk 
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There are two different regions of integration for k: 

(i) The UV region where all three fluctuations are short. This contributes at cubic orders in 
an expansion of the short-scale fluctuations. The dependence of these short scales on long- 
wavelength fluctuations is captured by a primordial trispectrum. Since we are working at first 
order in non-Gaussianity, we will not keep track of these contributions in the rest of this appendix. 

(ii) The region where one of the modes is a long-wavelength fluctuation. This happens either for 
|fe| < A or for \k — p — ^k\ < A. We assume that these two regions do not overlap substantially 
and therefore consider them as disjoint. (The error made in using this assumption is captured 
by higher-derivative terms). Since both of these regions give exactly the same contribution to 
cr'-l, we focus on fc < A and multiply the result by two. 

Extracting the long-wavelength mode, the first-order contribution to becomes 


4 



(1) 



d^fc 

(27r)3 


d^T) 

j^n^^p,k)F,{k,p + lk-k) 

X 6i{-p + lk)di{p + ^k- k) 


6i{k) . 


(B.9) 


Next, we shift the variable of integration, p p + ^k, and express the two short fluctuations i5i in 
terms of the initial potential ipsiTin)- Expanding the integrand in the limit of large p, we get the 
answer in (B.3). □ 


It is instructive to go back to position space. Eq. (B.3) then becomes 

(T*-^'(®,r) = cI\x,t) + cY{x,t)5{x,t) + kiix,T) d''d’-(p{x,T) d-, (B.IO) 

which is to be compared to the structure of the stress tensor in (3.5). Up to corrections of order 

we have 

= Ca[d^4)s{q)d^4>s{q)]A , (B.ii) 

c^^ki{x,T) = -^^[d^d^(ps{q)dkdid~‘^Mq)]A + ^ {[d'^s{q)dkMq)]ASj -b perms) , (B.12) 

where q{x,T) is the Lagrangian position of the Eulerian coordinate x, a € {0,1} and Ca are 
unimportant numerical factors. Let us make a few comments: 

• We see that the coefficients are in general non-local functions of (ps — in this example, this 

is the case only for the coefficient Since the inverse Laplacian in (B.12) acts only 

on short-scale fluctuations, ki is a non-local function of (ps within a region of size A“^ 
centered around the Lagrangian coordinate q. 

• Some terms such as 6 and seem to be missing in (B.IO). This is just an artefact of 
working to lowest order in perturbation theory, where these terms are indistinguishable 
from 6 and d^d^cp. 

Eq. (B.3) shows how the late-time a'^^{x,T) depends on the long-wavelength fluctuations and 
the initial short-scale fluctuations. So far, this is very general, in the sense that we haven’t 
specified the initial conditions for the short scales. Next, we will determine how these initial 
fluctuations are modulated by long-wavelength fluctuations in the presence of PNG. 
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B.2 Non-Gaussian Initial Conditions 


To determine the dependence of (Jp{k) on the long-wavelength fluctuations, we average (B.3) over 
the short scales. This boils down to replacing cjp(fc) by {ap{k))s. We then substitute (3.10) into 
the expression of {ap{k))s. The result receives contributions from different spins 


{ap{k))s = ap{k) + /nl [(cTp(fc))i°] + (CTp(fc))Pl -P (o-p(/c)) W] , (B.13) 

where ap{k) oc 5{k,) is the Gaussian contribution (which is independent of the long-wavelength 
fluctuations) and ((Tp(fc))i^^ represent the spin-L contribution 

(o-p(fc))[°] = {n/p)^P^{p) X V'(fc) , (B.14) 

(np(fc))P] = {p/p)^P^{p) (pipj) X V’*^(fc) , (B.15) 

(np(fc))W = ik'/p)^ P<t>ip) (PiPjPkPi) X , (B.16) 

where the fields V’; ^md were defined in eqs. (3.13), (3.21) and (3.25), respectively. As 


explained in the main text, higher-spin contributions are captured by higher-order tensors which 
can only contribute at higher order in the fluctuations and/or derivatives. For this reason, we 
did not consider them in this work. 


B.3 EFT Operators 

We split (T®-^ into Gaussian and non-Gaussian contributions 

= + (B.17) 


To obtain the Gaussian contribution we replace {ap{k))s by its Gaussian contribution ap{k). The 
non-Gaussian part gets contributions from the different spins: 




(B.18) 


To compute the spin-L contribution to (cr®-^)^, we replace cjp(fc) in (B.3) by (cTp(fc))i^^. Going to 
real space, we find 


(.T«)[»l = 1/3 


1 ^ 

^ + ^T9®di$ 


6 


7 / 


7' 


oo 

where we have defined the coefficient 


+ —¥^6 - —¥^dkdi^S^^ + 

49 49 49 


,3(a.a). 


/A 27r2 {p/p) 


A ■ 


(B.19) 

(B.20) 

(B.21) 

(B.22) 


We see that all terms consistent with the symmetries are generated in the effective theory — see 
eqs. (3.31), (3.32) and (3.33). 
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C Perturbation Theory and Connterterms 


In this appendix, we collect a few results from standard perturbation theory (SPT) and derive 
explicit expressions for the one-loop counterterms in the EFT-of-LSS. 


C.l Equations of Motion 

The equations of motion for the density contrast 6 and the velocity divergence 0 are 

drS + 0 = Sa , 

(dr + ?{)0 + S -\- Tg , 

where tq = — dj , and Sa^fi are nonlinear source terms 

Sa(k,T) = - / a{p,k-p)0{p,T)6{k-p,T) , 


, X ki- {ki + k2) 
a(ki,k2) = - 


Sf}(k, t) = - /3(p, k - p)0{p, T)0(k -p,t) , /3(ki, fca) = 

Jp 


kl 

(ki + ^ 2 )^ • k2 


2 

Using the scale factor a(T) as the evolution variable, the equations of motion become 

Vs 


(C.l) 

(C.2) 

(C.3) 

(C.4) 




-a^dl + ( - 3 ) ada + 




-\-a^d‘^ + ( 4 — -Vm ) ada + (2 — 311^) 


6 = Sp + Tg - nda(aSa) , 


0 = da{aS/3 + are) - -Q.m'HSa ■ 


(C.5) 

(C.6) 


Ve 


The linearized equation of motion for 5 (obtained by setting the right-hand side of (C.5) to zero) 
has the following growing mode solution 

5(i)(fc,a) = Di{a)6i{k) , (C.7) 

where 6i{k) describes the initial condition and Di(a) is the linear growth factor 


i( ) 2^°^-a H3(a0 • 


(C.8) 


In Einstein-de Sitter, the result in (C.8) reduces to Hi (a) = ajuin- 

Solving (C.5) with a delta-function source, doia — a'), gives the Green’s function for the 


evolution of 5: 


n ( cs( 2 1 Hi(a') 

Gs(a,a) = B(a — a )- 


D-(a) Hi (a) 


(C.9) 


o' [D.(a') Hi(a')J ’ 

and H_ = 'H/(a'Ho) « H^ A similar Green’s function Gg(a,a') exists for 0, but it won’t be 
needed in this work. To a remarkably good approximation, we have 

2 


do' Gs{a,a')n^{a')r{a')[Di{a')] 


iDi(a)r , 


(C.IO) 


(n - l){2n + 3)' 

where / = din Di/din a. The result in (C.IO) will be useful below. We have checked that this 
approximation is valid to better than 2% accuracy, for the values of n considered in this paper. 
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C.2 Perturbative Solution 

For 5(1) < 1, the solution to the nonlinear equations can be written as a series in powers of the 
initial dark matter contrast 5i (and integrals over the Green’s function Gs)- 

OO OO 

5(fc,a) = ^5(„)(fc,a) , 9{k,a) =. (C.ll) 

n=l n=l 

For a non-vanishing stress tensor in (C.2), the n-th order solution can be written as 

5(„)(fc,a) = 5SfT(fc,a) + 5[=„)(fc,a) + 5(^„)(fc,a) , (C.12) 

where 5®^"^ is the SPT result obtained with tq = + Tn = 0, 5^ is the solution sourced by the 

viscosity part of the stress tensor while 5'^ is sourced by the noise component of the stress 
tensor r^. We first review the SPT solution (§C.2.1) and then derive the expression for the 
one-loop counterterms (§C.2.2 and §C.2.3). 


C.2.1 SPT Solution 

The n-th order SPT solution can be written in terms of a convolution of the Green function (C.9) 
and lower-order SPT solutions. Using (C.IO), one finds that the n-th order solution is proportional 
to the n-th power of the linear growth factor T^i: 


Sf^)^{k,a)^D^{a)5n{k) . (C.13) 

A similar result holds for the velocity divergence 

df^^^{k,a)^DUa)dM ■ (C.14) 

The initial conditions 5n and On can be written as a convolution of powers of 5i: 

= f f {27rf3D{k-ki,„n)Fn{ku...,kn)di{ki)...6i{kn) , (C.15) 

J t/ kji 

Of^{k)= [ ... [ {27rf 6D{k - ki,„n) Gn{ki,..., kn) 6i{ki)... 5i{kn) , (C.16) 

J k\ J kfi 


where = ki + ■ ■ ■ + kn. Explicit expressions for the kernel functions Fn and Gn can be found 

in [2]. For instance, Ti = Gi = 1, while the (symmetrized) second-order kernel functions are 


F2ikr,k2) = ^ + l 
G2{kl,k2) = ^ + l 


fki-k2 ki-k2\ 

V kf + fei J 

fki-k2 ki-k2\ 

V kf + ) 


2{ki-k2f 
^ 7 klkl ’ 
4 (fci • k2f‘ 

7 k^k"^ 


(C.17) 

(C.18) 


C.2.2 Viscosity Counterterms 

In order to renormalize the one-loop bispectrum, we need to compute 5^ and 6^ up to second 
order. This requires knowing the stress tensor up to second order. First, let us focus on the 
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viscosity contribution to the stress tensor 


= -(f{a)A6^i) - 5r(a)/NL A^(i) , 

tP = -(f{a)A6^2) -ei(a) A((5 (i)) 2 - e2(a) A (s(i))^ - e3{a)di{s"^^^dj6^i)) 


(C.19) 


~ /n 


5'(a)[A^'(2) - 9i(5(i)9*^'(i))] +fi(i(a) A(^'(i)(5(i))+5r2(a)9i9j(sJ^^)^(i)) , (C.20) 


where 5 (i, 2 ) and refer to the SPT solutions and not the full solu¬ 

tion (C.12). We have defined the parameter d? = c^ + /(c^js + c^js) as the sum of the sound speed 
and the viscosity parameter (for more details, see [ 8 , 9]). Note that the velocity divergence 9 
does not appear in these equations because at second order it is completely degenerate with 5. 
The field 'I'(a;,r) = 'ijj(q{x,T)) admits an expansion in powers of the fluctuations, with = -0 
and T( 2 ) = 'V'lp ■ V^h; cf. eq. (4.1). 

The n-th order counterterms 6f , can be written as 

(n) 

6A{k,a)= / ... {27rf6D{k-ki,„n) F^{ki,...,kn\a)5(^i){ki,a) ...6(^i){kn,a) 

Jki Jkn 

+ /nl / ■■■[ {27rf6D{k-ki,,,n)H^{ki,...,kn\a)'il;{ki)...6(^i-){kn,a). (C.21) 

O J ^TL 

where and F[^ are kernel functions that are to be determined. In what follows, we will compute 
these kernels up to second order. 


First order 

The first-order counterterm is the solution to Vs and can therefore be written as 

^{i)ik,a) = -|(a)fe2(5(i)(fc,o) -7(a)/NL/c^4'(i)(fc) , (C.22) 

where we have defined 

= ~ n\ \ i Gs{a,a)(f{a')Di{a') , 

Di{a) Jai„ 

■y{a) = - f da Gsia,a') g{a') . 

J din 

We see that the parameters ^ and 7 will depend on the time dependence of d^{a) and g{a). For 
the one-loop power spectrum, this is not very important, since in the end we just need to fit the 
values of ^ and 7 at a given redshift. However, as we shall see in the next section, the second-order 
counterterm < 5 ^ 2 ) depends on the time dependence of these parameters. A convenient ansatz for 
the time dependence of the parameters in (C.23) and (C.24) is 

d\a) = [n{a)f{a)f[D 3 {a)r^+^ d^ , (C.25) 

g{a) = [H{a)f{a)f[D 3 {a)r^^^ g , (C.26) 


(C.23) 

(C.24) 
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where d and g are constants. Using (C.IO), we then have 

2 


?(o) = 




[md + l){2md + 7) 


(C.27) 

(C.28) 


Second order 

At second order, things are a bit more complicated, since the second-order solution will depend 
on the precise time dependence of the first-order solution Consider the equation of motion 
for 

Vs <5[=2) = Sf + ri^) - ndaiaS^^^) , (C.29) 

(2) (2) 

where Sa ^ and Sy ' are obtained by replacing one of the <5 or 0 in the convolution by their linear 
SPT solution and the other by the corresponding linear counterterm. Using (C.13) and (C.22) 
for and respectively, and replacing by and by the solution for 

(5^2) can be written as 


<5[^^(A:,a)= / F|(p, A: - p|a) 5(i)(p, a)(5(i)(fc - p, a) , 


=/nl / H 2 {p,k - p\a)ip{p)6(^i){k - p,a) . 


(C.30) 

(C.31) 


We have separated the solution into a Gaussian part and a non-Gaussian term We 


look at each of these terms in turn. 


d 2 ) 


• Gaussian contributions 

The kernel functions were computed in [8, 9, 36]: 

fcsja) = F2"(fci, fcsja) + Ff{ki,k2\a) + F^iki, fcsja) , (C.32) 

where F^ is sourced by the nonlinear terms in F!^^ is sourced by and F^ is 
sourced by the second-order SPT solution 5^2^"'" which appears in the stress tensor (C.20) 

as D — A (5^2^"'". 


o The kernel iTJ can be written as 

3 

F2{ki,k2\a) = - '^ei{a)Ei{ki,k2) , 
i=l 


where the time-dependent coefficients are 


ei(o) = - 


1 


[Di{a)\^ 


da' Gs{a, a') [Di{a')Y ej(a') 


(C.33) 


(C.34) 
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and the momentum kernels are 
Ei[ki, k2) = ki2 5 


E2{ki,k2) = k^ 


12 


(fci • k2 

'^ 1^2 


Esiki, k2) — 


1 2 1 

- k-\n —fcl * /b2 

6 12 -r 2 i 


ki2 ■ k2 ^ ki2 ■ ki 


kl 


kl 


o The kernel E!^^ can be written as 


F!^P{ki,k2\a) = -i{a)E^p{kuk2) , 


where ^(a) was defined in (C.23), and 
1 


^apiklik2') — 


2md + 9 

2mrf + 7 


+ 


2^{ki,k2){kl + kl) 

a(fci, ^ 2 ) (k2 + + 2)A:^) + {1 •H' 2}^ 


2{md + 2 ) 

o Finally, the kernel E 2 reads 

F2^(fci,fc2|a) = -i{a)Ed{ki,k2) , 

where we have defined 

E,(k. k,) - i^d + mmd + 7) 2 pj, , , 

Es{ku k2) - ^ 2)(2m, + 9) ’ 

and E 2 is the SPT kernel (C.17). 

Non-Gaussian contributions 

The kernel function can also be written as a sum of terms 

77|(A:i, k2\a) = TfJ(^i; ^ 2 |a) + H2^{ki, k2\a) + Ef^iki, k2\a) , 

where TfJ and have the same meanings as before, and H 2 is sourced by 
proportional to g{a) in the stress tensor. 

o The kernel can be written as 


H2{ki,k2\a) = - ^7i(a)Gi(fci,fc2) , 


i=l 


where the time-dependent coefficients are 

1 


7 i(«) = -- 


Diia) 

and the momentum kernels are 


da' Gs{a,a') Di{a')gi{a') , 


Gi(fci,fc 2 ) = 


12 ) 


G2{k^M) = 


{ki2 ■ k2)" 

kl 


--fc ?2 

3 


(C.35) 

(C.36) 

(C.37) 

(C.38) 

(C.39) 

(C.40) 

(C.41) 

(C.42) 
the term 

(C.43) 

(C.44) 

(C.45) 

(C.46) 
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(C.47) 


o The kernel can be written as 


{ki,k2\a) = --f{a)Gai3{ki,k2) , 


where 


Gaisiki, ^ 2 ) = 


2ma + 7 


+ 


2m„ + 5 


[{nig + l)a{ki, ^ 2 ) + a{k 2 , fci)] k\ 


(C.48) 


{nig + l){2mg + 7) 
o Finally, the kernel 77^ reads 

fc2|a) = -^{a)G^{ki, k2) 

where 


G\s,{ki,k2) = 


nig{2nig + 5) 


k 


2 ki • k2 


- ki2 ■ ki 


{nig + \){2nig + 7) 

C.2.3 Noise Counterterms 

Up to second order, the noise contributions to the stress tensor are 
= -didjJl^ , 

rP = di{6dj4^) - didj{Jl^6(^^)) - didj{J2^kisf[)) - . 


(C.49) 

(C.50) 


(C.51) 

(C.52) 


The re-th order noise counterterm 6’/^^ can be written as the sum of Gaussian and non-Gaussian 
contributions 


(n) 

It, 


S(^^{k,a) = 6-l;^^{k,a) + S(;^^^{k,a) . (G.53) 

The non-Gaussian part receives contributions from both noise terms in the initial conditions (3.8) 
and the noise terms in the stress tensor (3.35). We will derive explicit expressions for the coun¬ 
terterms up to second order. 

First order 

Solving Vs = Tn\ gives the hrst-order Gaussian solution 

{k,a) = kikj f da'Gs{a,a') Jq-'{ k,a) = kikjN'J {k,a) . (G.54) 

J CLln 

The non-Gaussian contribution to 5^^ only comes from the initial conditions. In fact, the solution 
is obtained by replacing <5(i)(fc, a) 1 —)■ f^i^M{k)il)j{k) in the SPT expansion. At first order, we get 


6-l[^^{k,a) = UM{k) Mk) 


(G.55) 
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Second order 

The second-order Gaussian contribution is obtained by solving (C.5) with replaced by the 

( 2 ,') 

Gaussian contribution to Tn and by replacing one of the 6 (or 0) in by their linear SPT 
solution and the other by the first-order Gaussian noise contribution (or ). Since we 
will not require the second-order Gaussian noise counterterm for this work, we will not explicitly 
compute it here and refer the reader to [8] for an explicit expression. 

The second-order non-Gaussian contribution is obtained in a similar way. The contribution 
coming from the initial noise term ?/)j is obtained by replacing e-)■ f^i^M{k)'ijjj{k) in the SPT 

expansion. Adding the term sourced by the non-Gaussian contribution in (G.52), we get 


S(2f^{k,a) = I kikjNll{k-p,a)'ip{p) + 2F2{k-p,p)5(i){k-p,a)M{p)'il}j{p) , (G.56) 


\ 2 ) 


where 




da' Gs{a,a')f^{k,a!) . 


(G.57) 
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D IR-Safe Integrands 


The one-loop integrals can have divergences as the loop momentum p or one of the external 
momenta ifcj approach 0. Although individual diagrams can be IR divergent, the equivalence 
principle guarantees that the sum will be finite. To avoid delicate cancellations of large integrals 
(which might affect the precision of the numerical computation) it is useful to dehne integrands 
which are well-behaved in the IR. For Gaussian initial conditions, this was done for the power 
spectrum in [7] and for the bispectrum in [8, 9]. In this appendix, we extend these results to 
non-Gaussian diagrams. 

D.l Power Spectrum 

Let us hrst consider the non-Gaussian contribution to the one-loop power spectrum (4.3) 

Pi2{k) = 2 f F2{k - p,p)Biii{k,\k - p\,p) = j pi2{p,k) . (D.l) 

J p J p 

The kernel function F 2 is divergent when either of its argument vanishes [2], cf. eq. (G.17). This 
means that the integrand is divergent in the limits p —)• 0 and p —)• fc. However, these IR 
divergences are unphysical and cancel in the integral. To make this cancelation manifest at the 
level of the integrand, we will hrst map the divergence at p = A: into a divergence at p = 0. To 
do so, let us split the region of integration as follows 

Pi2{k)= J Pi2{p,k) + I pi2{p,k) . (D.2) 

p<\k—p\ p>|fc—p| 

Changing the integration variable in the second integral, p 1 —)• fc — p, we get 

Fuik) = 2 I pi2(p, k) 0(|fc - p| - p) , (D.3) 

Jp 

where 0 is the Heaviside function. Although this integrand no longer has a divergence at p = k, 
there is still a divergence at p = 0. To remove this divergence, we hrst notice that 

hinF 2 (fc-p,p) = -hO(p°) . (D.4) 

P^O 2p^ 

This means that the problematic divergence disappears if we make the integrand invariant under 
the exchange p 1 —— p. We therefore write 

Pi2{p,k) = pi2{p,k)Q{\k - p\ -p) +pi2{-p,k)e{\k + p\ - p) . (D.5) 

This integrand is IR-safe: it is well behaved in the limits p — )• 0 and p — Fk. 

D.2 Bispectrum 

In the main text, we have written the non-Gaussian one-loop bispectrum as 

Aoop = [ [4ii(P> ^i) + ® + ^ 122 (P> ki) + 6 ( 22 (P, ki) + perms] , (D. 6 ) 

Jp 
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where 


b^^l^{p,ki) = ‘iF:i{ki+p,-p,k2)Biii{ki,p,\ki+p\)Pii{k2) , (D.7) 

4n(P:^i) = 3F3 (/ci,p, -p)Siii(/ci,A;2,fc3)Pii(p) , (D.8) 

422(p,^i) = 4F2(A:3 +P, -p)P2(p,fc2 -p)-Biii(^i, 1^3 +p|, 1^2 -p|)Pii(p) , (D.9) 

b^^ll{p,ki) = 2F2iki,k2)F2{ki - p,p)Biu{ki,\ki - p\,p)Pu{k2) ■ (D.IO) 


These integrands contain divergences when the loop momentum p approaches zero or ±fcj. Fol¬ 
lowing the same logic as for the power spectrum, we can write the bispectrum (D.6) as an integral 
of two IR-safe integrands 


= [ [b^'^\p,ki) + b^^\p,ki)] . (D.ll) 

Jp 

The integrands b^"^\p,ki) and b^^\p,ki) can be written as 

b^^\p,ki) = b['^\p,ki) + b^2^\p,ki) + bi^^\p,ki) , (D.12) 

b^^\p, ki) = b[^\p, ki) + b^2^\p, ki) , (D.13) 


where we have defined 

b['^\p,ki) = 6®2(^2 -p,fci) 0(1^2 -p| -p)0(|fe3 +p| - 1^2 -p|) 

-h 0(1^3 +p| -p)0(|fc2 -p| - 1^3 +p|) 
b 2 ^\p, ki) = 26f22(p> ki) Q{\ki - p\ - p) + 5 perms , 
b^^\p,ki) = 26^\\(p, fcj) Q{\ki +p\ -p) + 5 perms , 

b[^\p,ki) = ^ 6®2(P>^i)0(l^3 +P| -p)0(|fc2 -p| -p) + {p^ -p} 
b'^ 2 ^\p, ki) = bfll{p, ki) + 2 perms . 


-|- 5 perms , 


-|- 2 perms , 


(D.14) 

(D.15) 

(D.16) 

(D.17) 

(D.18) 


It is easy to check that the integrands do not have any divergences as p approaches one of 

the external momenta Fki. Moreover, while each of the individual integrands is divergent 

in the limit p —)■ 0, the sums and b^^'^ are finite. 
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E Notation and Conventions 


Symbol 

Relation 

Meaning 

Equation 

a 


scale factor 


T 

a dr = dt 

conformal time 


n 

= dln{a)/dT 

conformal Hubble parameter 


'Ho 


present value of % 


X 


comoving coordinate 


Q 


Lagrangian coordinate 


k 


long momentum 


P 


short momentum 




matter density in units of the critical density 




dark energy density 


h 


dimensionless Hubble constant 


P 


dark matter density 


6 

= bp/p 

dark matter density contrast 


0 

= diP 

velocity divergence 


<J(n) 


density contrast in SPT at order n 

(2.4) 

^(n) 


velocity divergence in SPT at order n 

(2.4) 

Fn 


kernel function in 

(C.15) 

Gn 


kernel function in 

(C.16) 

P 

^ mn 

= (<5(m)'5(n))' 

power spectrum in SPT 

(2.7) 

Blmn 


bispectrum in SPT 

(2.8) 

Ps 

= (M)' 

nonlinear dark matter power spectrum 

(2.7) 

Bs 

= {bbb)' 

nonlinear dark matter bispectrum 

(2.8) 



dimensionless power spectrum 

(5.23) 

Bs 


dimensionless bispectrum 

(5.24) 

</> 


Newtonian potential 

(2.3) 

$ 

A$ = b 

rescaled Newtonian potential 

(2.3) 

‘P 

(/ = T{k)Lp 

primordial potential 

(2.11) 

Pg 


Gaussian primordial potential 

(2.17) 

A<^(A;o) 

= 3.0 X 10-5 

amplitude of the primordial potential 

(2.13) 

ko 


pivot scale 

(2.13) 

T{k) 


transfer function 

(2.11) 

M{k) 


transfer function in the Poisson equation 

(2.12) 

Di 


linear growth factor 

(C.8) 

f 

= d In L)i / In o 

growth rate 

(3.31) 



primordial power spectrum 

(2.13) 

B^ 


primordial bispectrum 

(2.15) 

Us 


scalar spectral index 

(2.13) 

Sij 

= didj^ — ^bij A<1> 

tidal tensor 

(3.31) 
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Symbol 

Relation 

Meaning 

Equation 

V' 


correlation in the initial conditions 

(3.13) 



Eulerian definition of ip 

(3.17) 

ipj 


noise term in the initial conditions 

(3.8) 



spin-2 correlation in the initial conditions 

(3.21) 

^ijkl 


spin-4 correlation in the initial conditions 

(3.25) 



kernel function of the primordial bispectrum 

(2.17) 

ai 


coefficient in the Legendre expansion of Knl 

(2.20) 

A 


scaling dimension in 

(2.20) 

/nl 


amplitude of the primordial bispectrum 

(2.20) 

Pl 


Legendre polynomial of order L 

(2.20) 

Wa 


window function 

(3.1) 

Fa 


Fourier transform of Wa 

(3.1) 

rij 


stress tensor of the EFT 

(3.3) 

{P^s 


viscosity part of 

(3.7) 



noise part of 

(3.7) 

f3 


parameter in AP^ 

(3.35) 

Te 


EFT source in the Euler equation 

(4.20) 

Tv 


viscosity part of Tg 

(4.22) 

Tn 


noise part of Tg 

(4.23) 

Cs 


sound speed 

(2.10) 

^vis ? ^vis 


viscosity parameters 

(2.10) 

(f 

^ cl + ficlis + ctis) 

parameter in 

(4.22) 

^ii 9i 9i 


parameters in Tv 

(4.22) 



parameter in 

(4.31) 

7 


parameter in 

(4.32) 

e* 


parameter in 

(C.34) 

li 


parameter in 6 ^ 2 ) 

(4.40) 

Sa.p 


SPT quadratic source terms 

(C.3) 

^(n) 


viscosity counterterm at order n 

(4.26) 

7«) 


noise counterterm at order n 

(4.44) 

F^ 


kernel function in 6 ^^-^ 

(C.32) 

H^n 


kernel function in <5^^^ 

(C.42) 

Gs 


Green’s function for 6 

(C.9) 

Vs 


evolution operator in the fluid equation 

(C.5) 

Pli, 

= ('5(i)V')' 

correlation of (5(i) and pi 

(4.2) 

Pic 

- (^(i)<J(i))' 

correlation of (5(i) and 

(4.33) 

aHA) 


divergence in Pi 2 

(4.5) 



divergence in Pi 2 

(4.6) 
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Symbol 

Relation 

Meaning 

Equation 

a^{K) 


divergence in R 122 

(4.16) 



divergence in -B 122 

(4.17) 



velocity dispersion 

(4.14) 

pG 

^SPT 

= -Pii + Pi3 + P 22 

Gaussian SPT contributions to Ps 

(5.3) 

pNG 

-^SPT 

^Pl2 

non-Gaussian SPT contribution to Ps 

(5.4) 

pG 

^EFT 

= -2Cfc2Pn 

Gaussian EFT counterterm 

(5.3) 

pNG 

-'^EFT 

= -2'yk‘^Pi^ 

non-Gaussian EFT counterterm 

(5.4) 

PgPT 


Gaussian SPT contributions to Bs 

(5.5) 

nNG 

ZJgpT 


non-Gaussian SPT contributions to Bs 

(5.6) 

rG 

■Deft 


sum of Gaussian EFT counterterms 

(5.7) 

nNG 

-Deft 


sum of non-Gaussian EFT counterterms 

(5.8) 


= -B^pt + CB^ 

Gaussian SPT contributions plus B^ 

(5.14) 

Br 


non-Gaussian SPT contributions plus B^^ 

(5.14) 

B? 

= -Beft “ ^Bf 

sum of Gaussian counterterms {—B^) 

(5.16) 

BNG 

— rNG c rNG 

— -Deft 

sum of non-Gaussian counterterms {—B^^) 

(5.17) 
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